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INTRODUCTION 


Sea  ice  is  affected  by  ice  notion  in  two  wayB:  first,  the  velocity 
advects  ice  from  one  region  to  another,  and  second,  the  deformations  associ¬ 
ated  with  spatial  differences  in  velocity  are  important  to  the  heat  balance- 
through  open  water  forma tion--and  to  the  momentum  balance— through  ridging 
and  ice  stress.  This  report  describes  research  applicable  to  both  aspects 
of  ice  motion,  but  concentrating  on  the  spatial  structure  of  the  velocity 
field  and  what  this  structure  tells  us  about  deformation. 

The  work  is  described  in  three  parts.  The  first  is  the  measurement  of 
ice  displacement  from  synthetic  aperture  radar  (SAR)  on  the  SEASAT  satellite. 

The  important  feature  of  these  measurements  is  their  hf^h  spatial  density— 

2  km  Intervals  along  an  865  km  track— which,  clearly  shows  for  the  first  time 
the  motion  of  distinct,  rigid  pieces  of  ice  or  floes.  These  measurements 
show  the  difficulty — or  at  least  the  arbitrariness— in  defining  a  velocity 
derivative.  The  second  subject  is  the  analysis  of  these  data  and  FGGE  buoy 
data  to  provide  the  spatial  autocorrelation  function  of  the  velocity  field 
and  to  show  how  it  is  related  to  experiment  design  and  the  notion  of  (spatial) 
velocity  derivatives.  Last,  a  model  of  the  velocity  field  as  the  motion  of 
a  set  of  pieces  shows  the  uncertainty  in  estimating  the  amount  of  opening  or 
ridging  in  an  ice  cover  from  velocity  measurements  at  only  a  few  points,  as 
was  done  for  instance  during  AIDJEX. 

Two  papers  in  which  these  subjects  are  discussed  in  detail  are  presented 
as  appendices.  Appendix  A  is  a  paper  by  Hall  and  Rothrock  on  the  SAR  measure¬ 
ment  technique  and  data.  It  has  been  published  in  the  Journal  of  Geophysical 
Research.  Appendix  B  is  a  paper  by  Thorndike  on  ice  kenematics.  It  defines 
the  spatial  autocorrelation  function,  evaluates  it  from  the  SAR  and  buoy  data 
and  describes  the  model  of  the  piece-like  velocity  field;  it  will  appear  as 
a  chapter  in  a  NATO  volume  on  Air-Sea-Ice  Interaction. 

DISPLACEMENT  MEASUREMENTS  FROM  SEASAT  SAR 

The  objective  is  to  compare  two  images  of  the  same  ice  taken  at  two 
different  times  and  to  estimate  the  movement  of  the  ice  during  the  time 
interval.  The  difficulty  is  in  establishing  the  relative  location  of  the 
two  images.  Our  procedure  was  to  use  images  from  two  nearly  identical  orbits 
and  to  align  the  two  images  by  comparing  the  portions  over  land.  No  effort 
was  made  to  establish  absolute  position. 
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The  displacement  from  5  October  to  8  October  1978~or  equivalently, 
the  average  velocity  over  the  three  day  interval— was  measured  at  points 
•paced  about  every  two  kilometers  along  a  line  865  km  long  in  the  Beaufort 
Sea.  Figures  1  and  2  in  Appendix  A  show  the  location  of  the  images  and  the 
data.  Both  components  of  displacement  are  presented  as  a  function  of  dis¬ 
tance  along  track.  The  motion  of  individual  pieces  is  evident  in  the  figures 
and  can  in  fact  be  found  by  fitting  piece-wis';  rigid  body  velocity  profiles 
to  the  data.  The  movement  of  several  pieces  is  shown  in  Figure  3  of  Appen¬ 
dix  A. 

The  along  track  component  of  displacement  has  an  error  which  is  zero 
at  the  coast  growing  to  about  3  km  at  the  far  end  of  the  strip.  This  error 
is  due  to  stretching  in  the  optical  processor  and  could  be  reduced  by 
choosing  strips  with  land  in  both  ends  (our  strip  Included  Banks  Island  but 
stopped  short  of  Siberia).  The  across  track  component  of  displacement  has 
an  error  of  0.4  km.  This  error  could  be  reduced  by  use  of  orbital  data  and 
by  correcting  for  the  slant  range  distortion.  We  recommend  these  improve¬ 
ments  be  made  in  future  work  with  this  type  of  data. 

VELOCITY  FIELD  STRUCTURE 

The  autocorrelation  function  tells  hov  similar  the  velocities  are  at 
two  points  a  distance  r  apart.  If  r  is  as  small  as  a  few  kilometers,  the 
two  velocities  are  nearly  the  same;  so  their  correlation  is  unity.  At  large 
r,  the  two  velocities  are  unrelated— zero  correlation.  The  autocorrelation 
is  described  by  two  functions  of  r.  Suppose  the  two  points  lie  on  the  x-axis. 
B||(r)  is  the  correlation  of  the  x  velocity  components,  and  B|(r),  the  y 
components . 


B | | (r)  -  u(o)n(r) 


B^(r)  -  y(o)*v(r) 
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These  functions  have  been  evaluated  from  AIDJEX  and  FGGE  buoy  data  and  are 
shown  in  Figure  16  of  Appendix  B,  The  notable  feature  is  that  the  correla¬ 
tion  length  scale  is  about  1,000  km. 

These  functions  have  several  uses.  In  experiment  design  they  allow  one 
to  choose  an  economical  sample  spacing  precluding  either  redundant  data  or 
data  gaps.  They  can  be  used  to  estimate  interpolation  errors  for  a  given 
data  set  (see  15.3.3.2  of  Appendix  B).  One  can  estimate  from  them  statis¬ 
tics  of  deformational  quantities  involving  spatial  velocity  differences  (see 
§5.4.2  and  .3  of  Appendix  B). 

The  notion  of  the  velocity  gradient  du/dx  has  been  central  to  sea  ice 
modeling,  yet  there  is  some  difficulty  in  defining  this  quantity.  A  deriva¬ 
tive  is  normally  defined  as  a  limit  ,  but  the  SAR 

data  (Figure  12  of  Appendix  A)  show  that  such  a  limit  does  not  exist  every¬ 
where,  and  where  it  does  exist  it  describes  the  rotation  of  a  floe,  not  a 
"large  scale"  gradient.  The  derivative  exists  in  the  mean  squared  sense,  if 


11m 
h  o 


u(x  +  h)  - 


where  E  is  the  expected  value  operator, 
equivalent  to 


■  a  constant 

The  existence  of  the  limit  is 


variance  ^u(ac  +  h)  -  u(x)j  -  oha ,  as  h  o 

with  a  •  2.  The  SAR  data,  plotted  in  Figure  19  of  Appendix  B  give  a  *  1.3 
forw  (a  ■  1.1  for  v),  giving  further  evidence  that  the  velocity  is  too 
irregular  spatially  for  a  good  derivative  to  exist. 

The  lack  of  a  derivative  can  be  circumvented  in  two  ways.  The  average 
derivative  over  some  region  x  can  be  evaluated  as  an  integral  of  some 
velocity  component  around  the  boundary  of  the  region.  The  autocorrelation 
functions  then  allow  one  to  evaluate  the  error  introduced  by  estimating  that 
Integral  with  only  a  few  (buoy)  data  points.  (See  §5.4.5  of  Appendix  B.) 

It  is  important,  though,  to  remember  that  the  value  of  the  average  deriva¬ 
tive  depends  on  the  size  of  the  sampling  region. 
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Alternatively,  one  can  replace  the  notion  of  derivative  with  velocity 
differences.  (See  15.4.2  of  Appendix  B.)  Velocity  differencea  exist  with¬ 
out  question;  they  can  naturally  be  associated  with  strain  rates,  ar.d  their 
covariances  can  be  estimated  from  5| |  and  B^. 

MODEL  OF  THE  PIECE-LIKE  VELOCITY  FIELD 

The  velocity  field  is  demonatrably  the  movement  of  rigid  pieces.  A 
model  incorporating  this  feature  can  be  used  to  examine  the  practice  of 
estimating  opening  and  closing  ^ridging)  from  velocities  measured  at  three 
or  four  points.  We  take  the  pieces  to  be  defined  by  a  random  set  of  lines 
(a  Poisson  field);  the  average  piece  size  of  16  km  is  taken  from  the  SAR 
data.  The  differential  movement  of  the  pieces  is  given  a  Gaussian  distri¬ 
bution.  Many  realizations  of  this  motion  are  simulated.  For  each,  the 
total  opening  (and  closing)  at  the  floe  boundaries  is  plotted  versus  a 
"large  scale"  estimate  of  deformation  made  from  the  velocity  at  three  points. 
Each  realization  gives  one  data  point.  The  ensemble  of  points  is  shown  in 
Figure  24  of  Appendix  B.  The  scatter  shows  how  poorly  a  three  point  strain 
rate  estimate  allows  us  to  determine  the  actual  small  scale  opening  and 
closing.  A  better  parameterization  would  account  for  the  scatter. 

RECOMMENDATIONS  FOR  FUTURE  WORK 

Considerably  mo,'e  can  be  done  with  SEASAT  SAR. 

1)  Further  assessment  of  errors  and  tests  of  algorithms  for  eliminating 
errors  should  be  made  by  examining  pairs  of  strips  over  land.  Strips  a 
thousand  kilometers  long  over  the  Canadian  archipelago  would  be  suitable. 

One  would  treat  the  strips  as  if  they  were  over  ice  with  land  in  each  end, 
and  apply  all  corrections  for  slant  range  distortion,  orbital  differences 
and  stretching  by  the  optical  processor  aud  then  test  the  central  portion 
of  the  strip  to  assess  the  remaining  positioning  errors. 

2)  SAR  ice  displacement  measurements  should  be  made  to  extend  the  present 
measurement  back  into  the  summer  season.  One  might  expect  substantial 
seasonal  variations  in  the  sizes  of  pieces  and  the  variance  of  the  motion. 

3)  A  sequence  of  three  or  more  SAR  strips  should  be  analyzed  to  assess 
the  persistence  of  floe  boundaries.  Do  the  velocity  jumps  occur  at  the  same 
floe  boundaries  for  weeks,  or  do  floe  boundaries  change  after  several  days? 
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The  answer  la  relevant  to  paramaterising  opening  and  closing  in  terns  of 
sparse  velocity  aeasurements. 

Other  new  work  could  include: 

4)  The  spatial  velocity  autocorrelations  should  be  computed  for  the  1980 
Arctic  buoy  data. 

5)  The  randon  piece-like  velocity  field  model  predicts  considerable 
uncertainty  in  our  estimates  of  open  water  formation  and  ridging  based  on 
velocity  measurements.  The  implication  is  that  the  thickness  distribution 
and  heat  flux  estimates  have  a  similar  uncertainty,  which  has  not  yet  been 
evaluated. 
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APPENDIX  A 

SEA  ICE  DISPLACEMENT  FROM 
SEASAT  SYNTHETIC  APERTURE  RADAR 

by  R.  T.  Hall  and  D.  A.  Rothrock 

This  appendix  appeared  in 
the  Journal  of  Geophyaioal  Research % 
Volume  86(C11) ,  11,078-82 


ABSTRACT 


A- 1 


Images  obtained  by  a  synthetic  aperture  radar  on  SEASAT  have  been  used 
to  measure  sea  ice  displacements  over  a  three  day  interval  in  October  1978. 

The  position  of  a  natural  ice  feature  was  measured  on  5  October  and  again  on 
8  October;  the  difference  in  these  two  positions  is  a  displacement.  The  dis¬ 
placement  of  many  features  was  measured.  The  tracked  features  lie  roughly 
along  a  line  and  are  quite  dense— about  2  km  apart— over  a  distance  of  865  km. 
The  displacements  are  about  twenty  kilometers.  Displacement  errors  grow  with 
distance  from  shore  becoming  as  large  as  3  km.  The  graph  of  displacement 
versus  distance  has  occasional  discontinuities  of  several  kilometers.  Displace¬ 
ment  discontinuities  are  accurate  to  ±  0.07  km  along  track  and  3X  of  their 
magnitude  across  track. 

INTRODUCTION 

The  motion  of  sea  ice  varies  on  all  scales  up  to  10s  km.  The  larger  scale 
motions  are  driven  by  atmospheric  pressure  systems.  On  scales  less  than  about 
102  km  the  motion  is  affected  by  individual  ice  floes.  The  motion  has  tradi¬ 
tionally  been  observed  by  tracking  particles  of  ice:  a  couple  of  drifting 
manned  camps,  or  a  few  buoys  or  radar  reflectors.  Each  tracer  is  costly, 
limiting  any  experiment  to  on  the  order  of  ten  points.  Such  methods  do  not 
reveal  the  spatial  structure  of  the  field  of  motion.  What  is  needed  are 
records  thousands  of  kilometers  long,  sampled  every  kilometer  or  so. 

Tracking  natural  ice  features  in  satellite  imagery  improves  the  sampling 
rate.  Following  recognizable  leads  in  LANDSAT  visual  photographs  can  provide 
records  over  500  km  long  with  about  25  irregularly  spaced  data  points — an 
average  spacing  of  20  km  [ Nye ,  1975].  Synthetic  aperture  radar  (SAR)  has  the 
potential  to  provide  yet  more  dense  measurements.  Its  high  resolution  and 
sensitivity  to  surface  roughness  allow  individual  roughness  features  to  be 
tracked,  and  these  are  more  abundant  than  new  leads. 

The  last  decade  has  seen  a  considerable  amount  of  radar  imagery  collected 
from  aircraft.  Airborne  imaging  radar  data  from  1975  was  used  to  provide  ice 
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displacements  at  spaclnga  of  about  10  km  over  an  area  roughly  10*  km1 
[Bryan  at  al.  1977;  L*b«rl  at  al.  1979].  But  In  general  the  distortion 
caused  by  the  vide  look  angle  and  the  aircraft  motion  discouraged  quanti¬ 
tative  analysis.  It  was  not  until  generous  amounts  of  quasi-corracted 
Imagery  from  the  L-band  SAR  on  SEASAT  became  available  through  NOAA  that 
the  general  sea  ice  community  could  experience  firsthand  the  tantalizing 
potential  of  radar  data. 

This  paper  describes  a  technique  for  measuring  ice  displacements  from 
these  SAR  images.  It  is  not  intended  to  be  a  comprehensive  analysis  of  the 
performance  of  the  system,  but  rather  a  report  on  some  encouraging  results. 
Measurements  are  presented  of  both  components  of  displacement  versus  a  single 
horizontal  coordinate.  This  data  record  is  865  km  long  and  has  417  data 
points — roughly  one  point  every  2  km.  Measurement  errors  are  assessed  by 
applying  the  same  method  to  images  of  land  where  there  are  no  displacements. 
The  errors  are  discussed  in  terms  of  the  satellite  system  and  the  optical 
data  processing  procedure.  It  appears  that  SEASAT  SAR  can  provide  3  day 
displacements  accurate  to  several  kilometers  and  displacement  discontinuities 
accurate  to  one  or  two  hundred  meters.  With  sufficient  care,  these  errors 
could  be  reduced. 


THE  DATA 

SEASAT  operated  from  early  July,  until  10  October  1978,  collecting  SAR 
data  over  swaths  100  km  wide.  More  than  10s  km  of  arctic  data  were  collected, 
half  of  this  over  sea  ice  in  the  Beaufort,  Chuckchl  and  East  Siberian  Seas 
south  of  75 ®N.  Often  the  sensor  was  turned  on  over  North  America,  crossed 
the  Arctic  Ocean  and  was  turned  off  over  Siberia.  Hence,  thete  is  often  land 
at  both  ends  of  the  swath  which  can  be  used  both  to  eliminate  some  errors  in 
the  data  and  to  assess  those  errors  remaining. 
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The  novel  feature  of  SAR  la  that  It  coobinea  data  taken  at  different 
positions  along  the  flight  path  to  synthesiaa  a  radar  image  from  an  antenna 
much  larger  than  actually  used— thus  providing  the  high  resolution.  [See, 
for  example,  Harg*rt  1970,  or  Assvea,  1975.]  The  data  are  transmitted  from 
the  satellite  and  stored  in  digital  format.  The  synthesis  of  the  bulk  of 
this  high  resolution  data  is  presently  pet  formed  by  the  Jet  Propulsion  Labora¬ 
tory  in  Pasadena  by  an  optical  rather  than  digital  process.  The  100  km  swath 
is  processed  as  four  adjacent  25  km  swaths.  The  scale  is  1:500,000;  hence, 
two  millimeters  on  the  film  strip  represent  a  distance  of  one  kilometer.  The 
resolution  of  our  images  is  about  0.035  km  on  the  ground;  0.07  tmn  on  the  film. 

The  film  strips  are  70  mm  positive  transparencies.  Not  all  are  of  good  quality 
we  received  one  strip  with  a  systematic  blurring  and  elongation  of  features 
and  another  with  double  images  for  each  pixel.  The  technical  details  of  the 
SEASAT  SAR  system  are  reviewed  by  Jordan  [1980]. 

These  images  of  sea  ice  look  much  like  visual  photographs.  Many  familiar 
shapes  and  structures  can  be  seen:  leads,  floes,  and  ridges.  The  same  shapes 
are  clearly  recognizable  in  views  of  the  same  ice  taken  at  different  times. 

The  correspondence,  however,  between  radar  and  visual  images  is  not  always 
consistent  [see,  for  example,  bsyan,  1976]. 

To  measure  displacement,  one  needs  to  know  the  position  of  a  recognizable 
ice  feature  imaged  from  two  successive  orbits.  Since  there  was  no  control  on 
absolute  location,  we  chose  to  use  orbits  three  days  apart  (or  a  multiple  of 
three)  because  these  orbits  are  almost  identical.  We  examined  in  detail  orbits 
1396,  1439,  and  1482  on  October  2,  5  and  8,  knowing  from  NOAA  VHRR  images  that 
there  was  ice  movement  during  this  period,  and  wanting  the  most  vinter-like 
ice  conditions.  The  data  reported  here  are  drawn  from  the  latter  two  orbits. 

The  ground  swath  is  shown  in  Figure  1.  Furthermore,  we  looked  only  at  features 
roughly  in  the  center  of  a  single  25  km  swath,  thus  avoiding  any  need  to  register 
the  adjacent  film  strips  and  also  minimizing  the  effect  of  slant-range  distortion. 


MEASUREMENT  TECHNIQUE 


To  identify  features  in  two  successive  images,  and  to  find  their 
initial  and  subsequent  positions  in  a  common  coordinate  system,  we  have  used 
a  Bendix  Analytical  Plotter  Model  C.  This  machine  normally  functions  as  a 
stereographic  mapping  tool;  we  have  used  it  as  an  x,y  digitizer  with  the 
crucial  feature  that  two  images  can  be  viewed  simultaneously.  A  few  guide 
lines  are  drawn  on  the  images  to  help  keep  the  viewer  oriented  as  he  searches 
for  features.  The  two  images  are  placed  on  the  two  viewing  stages  of  the 
plotter,  and  the  stages  are  independently  moved  until  the  two  cursors  are 
located  on  the  same  feature  in  both  images.  The  x,y  coordinates  of  the  cursor 
positions  relative  to  the  viewing  stages  are  then  recorded.  With  this  equip¬ 
ment;  one  can  measure  about  twenty-five  data  points  in  an  hour. 

Since  the  Bendix  plotter  is  designed  to  accept  24  cm  x  24  cm  aerial 
photographs,  our  long  film  strips  had  to  be  analyzed  piecemeal  in  overlapping 
segments  representing  roughly  100  km  of  ground  coverage.  To  reassemble  the 
data  from  the  segments  into  one  long  data  record,  two  reference  points  were 
marked  in  the  overlapped  regions  of  each  pair  of  adjacent  segments  and  measured 
in  both  reference  frames.  Thus,  the  data  from  each  segment  could  be  trans¬ 
lated  into  the  reference  frame  of  the  previous  segment,  and  the  segments 
reassembled  into  one  long  record. 

The  coordinate  system  is  defined  by  a  combined  use  of  time  marks  and  land. 
The  time  marks  are  dots  spaced  along  the  edges  of  the  film  strip.  They 
define  both  the  inboard  and  outboard  edge  of  the  25  km  ground  swath  even  though 
the  imagery  may  extend  a  bit  outside  them.  They  would  also  define  distance 
along  track  except  for  some  glitch  in  the  system.  Our  x-axis  is  taken  to 
run  through  the  inboard  time  marks  on  both  images.  The  y-axis  is  perpendicular 
to  x  through  the  most  westerly  land  feature  on  Banks  Island. 


A-' 6 


Ice  displacements  (u,v)  from  initial  positions  (x,y)  are  plotted  in 
Figure  2  versus  x  only.  The  data  disclose  the  motion  of  Individual  pieces 
of  ice.  Each  piece  has  a  solid  body  rotation  and  translation;  between 
pieces  are  jumps  in  displacement.  Figure  3  illustrates  this  by  showing  a 
line  drawn  from  the  first  to  the  last  point  on  each  piece  both  before  and 
after  the  displacement.  We  see  rotations  of  several  degrees,  translations 
of  ten  to  twenty-five  kilometers,  and  displacement  jumps  up  to  several 
kilometers. 

These  data  resolve  the  discontinuous  nature  of  the  field  of  motion,  and 
they  do  so  over  a  great  enough  distance  to  show  many  jumps.  Thus  we  now 
have  a  data  set  from  which  to  theorize  about  the  structure  of  the  ice  dis¬ 
placement  field  and  how  the  discontinuities  in  it  change  the  ice  cover.  The 
opening  and  closing  of  leads,  however,  cannot  be  obtained  directly  from  this 
one  dimensional  slice  of  a  two  dimensional  displacement  field.  They  require 
the  component  of  a  displacement  jump  normal  to  a  line  of  discontinuity  (or 
crack  or  lead),  and  we  have  not  measured  the  directions  of  these  lines.  In 
fact,  these  directions  are  not  well  defined  on  the  images. 

Plots  of  the  quantities  shown  in  Figure  2  were  made  to  the  same  scale 
as  the  image  to  see  what  correspondence  there  was  between  leads  and  displace¬ 
ment  jumps.  The  leads  of  this  October  scene  are  not  the  meandering  angular 
lines  of  April  ice,  but  are  more  lacy  like  the  open  water  surrounding  summer 
floes.  These  lacy  leads  change  shape  and  size  with  deformation.  One  might 
have  expected  the  big  events  of  deformation  to  occur  where  there  was  some 
striking  lead-like  feature,  but  this  was  generally  not  the  case.  At  all 
large  deformations  we  could  see  changes  in  leads,  but  the  leads  were  usually 
modest  (less  than  a  kilometer)--not  those  that  stood  out  for  their  size  (2  to 
A  km).  Some  areas  which  had  many  leads  and  looked  ripe  for  deformation  moved 


We  tried  to  find  a  measurable  feature  at  least  every  2  km  in  x.  In 
some  relatively  featureless  stretches,  this  was  difficult;  in  others,  we 
could  have  found  many  data  points  in  2  km.  Often  one  chooses  a  feature  in 
one  image  but  has  difficulty  finding  its  counterpart  in  the  other  image. 

There  are  two  effects.  In  successive  orbits,  a  moving  object  is  viewed  from 
different  angles  from  which  it  has  a  different  appearance  to  the  radar.  In 
addition,  the  ice  features  in  question  are  at  the  limit  of  the  resolution 
and  are  thus  defined  by  a  different  grouping  of  pixels  (picture  elements)  in 
successive  images.  Hence,  a  small  feature  prominent  in  one  image  may  be  absent 
in  the  other  or  distorted  enough  to  make  one  uncertain  where  the  cursor  should 
be  placed.  A  somewhat  lower  magnification  than  the  six  power  eyepiece  in 
our  Bendix  plotter  would  actually  help  one  Identify  features.  On  land, 
features  were  relatively  easy  to  identify  and  measure,  presumably  because  the 
land  is  stationary  and  we  had  chosen  to  look  at  nearly  identical  orbits. 

With  the  cursor  on  an  ice  feature  in  one  image,  our  uncertainty  in  locating 
the  feature  in  the  other  image  is  comparable  to  the  resolution  ±  35  m.  This 
error  pertains  to  positioning  a  bright  return  in  the  film  image.  There  is 
still  the  question  of  how  closely  a  bright  feature  seen  in  two  images  corres¬ 
ponds  to  one  particular  ice  feature.  We  have  not  tested  this,  but  the 
positioning  error  for  ice  features  can  only  be  greater  than  that  for  image 
features. 


OBSERVED  DISPLACEMENTS 

The  25  km  wide  strips  for  the  orbits  we  examined  (Figure  1)  include  460 
km  over  the  Canadian  archipelago  (x  <  0)  and  865  km  of  drifting  ice  between 
x  ■  285  and  x  *  1150  km.  We  made  no  measurements  in  the  near  shore  region 
(0  <  x  <  285)  where  large  displacements  carried  the  ice  to  the  extreme  edge 
or  completely  out  of  the  second  image. 
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rigidly.  Areas  which  moved  approximately  rigidly  but  did  show  some  relative 
displacement  (several  hundred  meters) — like  the  noisy  segment  between  x  ■  385 
and  x  *  435  km — do  in  fact  show  many  sites  at  which  this  deformation  could 
have  occurred.  This  sort  of  motion  seems  consistent  with  an  autumn  pack 
with  many  small  floes  being  consolidated  into  a  winter  pack  with  larger  pieces. 


ERRORS 

Our  interest  is  in  the  accuracy  of  relative  position  (x»y),  of  displacement 
(u,v)  and  of  displacement  differences  (Au,Av).  We  have 


u  *  x2  -  x1 
a  a 


2  1 
v  ■  y  -  y 

a 


(1) 


Au  *  (x2  -  x*)  -*  (x2  -  x1) 
b  b  a  a 

Av  -  (y2  -  y>)  -  (y2  -  y*) 

where  the  subscript  identifies  the  ice  particle  and  the  superscript,  a  time. 
Discontinuities  in  displacement  occur  for  only  a  few  unique  pairs  of  adjacent 
particles. 

Systematic  errors  arise  from  various  sources  which  we  will  discuss  in 
turn  below.  We  cannot  now  eliminate  these  errors  from  tha  data;  the  discussion 
is  aimed  at  eliminating  them  from  future  measurements  of  these  SAR  images. 

There  is  a  random  zero  mean  error  due  to  the  35  m  resolution  of  the  system. 

Its  standard  deviation  is  35  m  in  (x,y),  /T" times  this  or  50  m  in  (u,v)  and 

70  m  in  (Au,Av).  It  is  seen  in  Figure  4  as  the  vertical  scatter  about  the  trends. 


Along  track  stretching.  The  mechanics  of  the  processor  can  introduce  stretching 
of  up  to  IX  in  the  along  track  direction.  We  have  compared  our  Images  to  maps 
and  find  0.25  to  0.50*  stretching.  Conservatively  then,  we  take  our  values 
of  x  to  have  an  error  of  0.005  x. 
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When  two  Images  are  compared,  this  stretching  may  add  or  partially  cancel. 
We  find  that  the  image  on  October  8th  is  compressed  about  0.26%  relative  to 
that  on  the  5th.  This  is  the  trend  in  Figure  4.  Since  this  trend  is  not 
necessarily  maintained  out  over  the  ice,  we  do  not  correct  the  data.  If 
we  had  fixed  land  points  in  the  Siberian  end  of  our  images,  we  would  probably 
have  enough  information  to  remove  or  reduce  this  error.  The  resulting  error 
in  u  is  roughly  -0.0026  x;  at  the  far  end  of  the  strip  (x  ■  1150  km)  this 
amounts  to  3  km.  For  displacement  differences  Au,  the  errors  in  neighboring 
u’s  effectively  cancel. 

Across  track  drift.  Figure  4  shows  a  nearly  constant  error  in  v  of  -0.0001  x 
+0.38  km.  As  with  u.  this  error  cancels  when  one  calculates  Av  for  neighboring 
features.  The  error  in  v  can  arise  from  several  effects,  two  related  to  the 
orbits  and  the  third,  to  the  processor.  The  two  orbits  are  not  quite  identi¬ 
cal  but  cross  at  a  slight  angle  near  their  tops.  This  causes  an  apparent 
across  track  displacement  which  grows  with  distance  from  the  crossing.  The 
orbital  data  needed  to  remove  this  error  are  available.  The  other  orbital 
contribution  to  this  error  is  made  by  the  practice  of  taking  satellite  alti¬ 
tude  to  be  fixed  when  it  actually  is  slowly  changing  along  the  orbit.  Enough 
is  known  about  the  actual  altitude  to  remove  these  errors.  Occasionally 
(perhaps  once  or  twice  per  orbit)  the  assumed  altitude  is  corrected,  and  at 
these  locations  along  the  orbit,  there  is  an  apparent  jump  in  the  across  track 
position.  There  were  no  altitude  corrections  in  the  images  we  studied. 
Different  processor  set-ups  can  introduce  an  offset  in  the  across  track  loca¬ 
tion  of  the  orbit  on  the  order  of  a  few  hundred  meters.  This  can  only  be 
detected  and  removed  by  locating  the  x-axis  with  a  land  point  as  was  done 
with  the  y-axis.  Again  land  points  at  each  end  of  the  pass  (in  the  absence 
of  altitude  corrections)  would  remove  these  sources  of  error. 
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Across  track  distortion.  The  relation  between  the  slant  range  (the  distance 
from  the  satellite  to  a  point  x,y  on  the  ground)  and  the  across  track  coordi¬ 
nate  y  is  non-linear.  For  convenience  it  is  approximated  to  be  linear, 
introducing  a  distortion.  The  specification  for  the  optically  processed 
data  is  that  this  distortion  be  no  more  than  3%  of  the  distance  from  the 
centerline  of  a  25  km  strip.  We  have  no  independent  check  and  so,  use  this 
value.  The  greatest  distance  from  any  of  our  data  points  to  the  centerline 
is  6.6  km,  with  an  error  of  0.2  km.  The  3%  systematic  error  in  y  causes 
a  displacement  error  of  3%  of  v,  and  the  error  in  Av  is  3%  of  Av.  It  would 
be  a  simple  matter  to  correct  all  y  values  for  the  across  track  error. 

DISCUSSION 

Despite  some  early  pessimism  about  the  quality  of  the  SEASAT  SAR  optical 
data,  from  them  one  can  obtain  high  quality  measurements  of  ice  movement.  The 
major  errors  in  our  measurements  can  be  removed  with  more  work  and  information. 
One  needs  land  in  both  ends  of  the  ice  imagery,  orbital  parameters  (equatorial 
crossings),  altitude  data,  and  the  correct  relation  between  the  across  track 
coordinate  and  the  slant  range. 

There  is  no  handicap  in  working  with  images  which  have  not  been  digitally 
corrected.  Coordinates  of  features  can  be  measured  from  the  slightly  distorted 
Images  and  can  be  corrected  by  the  same  algorithms  one  would  use  to  correct 
the  whole  image,  with  considerable  saving  in  processing  cost  and  effort. 

This  statement  may  be  as  true  for  pattern  recognition  and  area  measurements 
as  it  is  for  displacement  measurements.  Others  have  made  this  point  [e.g. 
Leberlt  et  al.t  1979],  but  we  repeat  it  here  for  emphasis. 

There  are  many  ways  t.o  look  at  SAR  data.  Our  approach  has  been  limited 
to  making  observations  along  a  line.  Thus,  we  have  avoided  making  mosaics 
and  have  minimized  the  problem  of  mapping  different  orbits  into  a  common 
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coordinate  system.  Neither  hwe  we  attempted  to  establish  absolute  position 
except  for  illustrating  the  orbit  In  Figure  1.  Other  researchers  attempting 
a  more  general  treatment  will  encounter  more  difficulties  than  we  have.  We 
are  convinced  though  of  the  usefulness  of  these  data  and  the  pity  of  SEASAT's 
early  failure. 

There  are  more  pairs  of  congruent  strips  of  SAR  data  between  July  and 
October  1978.  Their  analysis  could  be  expected  to  show  some  dramatic  seasonal 
changes  in  the  velocity  field  structure,  both  in  the  sizes  of  rigid  pieces  and 
in  the  variance  of  the  motion.  There  does  not  seem  to  be  enough  data  to 
explore  regional  differences. 

Given  that  knowledge  of  ice  motion  is  essential  to  any  study  of  mass 
balance,  high  resolution  radars  are  a  pronising  source  of  data  on  an  opera¬ 
tional  level.  Their  disadvantage  is  their  high  data  rate.  An  alternative 
worth  studying  is  collecting  radar  data  only  in  a  small  fraction  of  the  polar 
oceans  in  data  windows  say  30  km  in  diameter  spaced  on  a  500  km  grid.  The 
data  set  would  serve  the  same  purpose  as  an  array  of  buoys.  The  details  of 
the  motion  would  not  be  observed  operationally,  but  their  statistics  would 
be  known  from  a  few  records  of  the  sort  reported  here. 
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October  1978.  Solid  lines  outline  the  portion  of  the  path  where 
displacement  measurements  were  m'.de. 

Figure  2.  Sea  ice  displacements  between  5  and  8  October  1978: 

(upper)  along  track  displacament  vs.  along  track  position) 

(lower)  across  track  displacement  vs.  along  track  position. 

Figure  3.  Motion  of  ice  pieces  from  5  to  8  October  1978.  The  line  segments 
represent  floes  moving  rigidly  from  the  lower  position  tc  the 
upper  position.  For  instance,  the  segment  AB  on  5  October  moved 
to  A*Br  by  8  October. 

Figure  4.  Apparent  displacement  over  land  from  5  to  8  October  1978.  Since 
the  true  displacement  is  zero,  these  data  represent  measurement 
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KINEMATICS  OF  SEA  ICE 


Kinematics  is  that  branch  of  physical  science  dealing  with  motion 
itself  and  not  with  the  causes  of  the  motion.  Attention  is  focused,  in 
this  chapter,  on  how  sea  ice  moves.  The  deeper  question  of  why  sea  ice 
moves  as  it  does  is  addressed  in  the  following  chapter  on  sea  ice  dynamics. 

For  some  scientific  and  practical  applications,  knowledge  of  the  ice  kine¬ 
matics  (as  opposed  to  the  dynamics)  is  sufficient.  As  examples,  consider 
gouging  of  the  sea  floor  in  shallow  water,  the  transport  of  natural  debris 
or  pollutants,  the  impact  loads  on  structures,  navigability  in  icy  regions, 
and  the  movement  of  the  ice  edge.  Aside  from  these  applications,  the  study 
of  kinematics  is  a  basic  step  toward  understanding  the  dynamics. 

Figure  1  illustrates  the  idea  that  the  motion  of  sea  ice  is  the  response 
of  the  ice  pack  to  external  forces.  To  understand  the  response  of  the  ice 
one  needs  to  monitor  both  the  driving  forces  and  the  notion.  Since  we  will 
discuss  only  the  motion  here,  there  will  be  the  underlying  ambiguity  of 
whether  the  results  inform  about  properties  of  the  ice  itself  or  about  the 
driving  forces.  We  will  see  for  example  that  the  ice  motion  is  nearly  non- 
divergent.  But  we  should  not  conclude  on  this  basis  alone  that  the  ice  is 
an  incompressible  material.  The  correct  explanation  could  be  that  the  npt 
driving  force  is  itself  nearly  non-dive rgent.  Because  of  this  ambiguity 
in  the  interpretation  of  kinematic  data,  it  is  useful  to  think  of  the  study 
'j  of  the  kinematics  as  part  and  parcel  of  the  broader  study  of  sea  ice  dynamics. 

The  fundamental  kinematic  notion  is  that  all  pieces  of  ice  have  an 
identity  which  is  preserved  in  time.  If  we  identify  a  piece  of  ice  say  by 
making  a  small  mark  on  the  ice  surface  at  position  X  at  time  o  ,  the  assumption 
is  that  at  some  later  time  -fe  ,  the  piece  of  ice  as  identified  by  that  mark  will 
be  at  a  position  X  .  This  defines  the  position  function  *(+,X)  with  . 

We  do  not  mean  to  imply  by  this  that  ice  floes  retain  their  identity  indefi¬ 
nitely,  only  that  it  is  possible  to  track  individual  points.  Two  points 
originating  on  the  same  ice  floe  may  very  well  wander  apart  in  time,  but  we 

assume  we  could  keep  track  of  them. 

There  are  difficulties  with  the  notion  of  a  position  function.  It  is 
natural  to  consider  the  space  domain  to  be  the  two  dimensional  surface  of 
the  Arctic  Ocean,  sayX*  At  4*o  there  will  generally  be  some  points  in  2 
which  are  not  covered  by  ice,  so  X^X)  is  not  defined  for  some  Xf2. 
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Further,  ice  is  always  being  created  and  destroyed  thermodynamically,  sc 
the  position  function  *(4,X,)  is  only  defined  for  an  interval  of  time. 

Finally,  the  idea  that  the  small  mark  retains  its  identity  is  question¬ 
able.  This  is  a  common  problem  in  the  definition  of  the  motion  of  a 
turbulent  fluid.  After  a  time  the  marked  fluid  has  mixed  with  th;  unmarked 
and  it  is  no  longer  possible  to  assign  a  position  to  the  mark.  Pomething 
of  the  same  sort  happens  with  sea  ice.  Despite  these  shortcomings,  the 
notion  of  a  position  function  is  the  best  we  have.  Furthermore  this  notion 
corresponds  exactly  with  most  of  the  observations  we  have.  We  will  use 
the  term  trajectory  for  the  function  *(4,X,)  where  X  is  held  constant.  The 
initial  position  X  serves  to  label  ice  particles  and  we  refer  to  the 
trajectory  x(4,X„)  of  the  particle 


5.1  Observations 

The  most  common  observations  of  ice  motion  are  of  trajectories  of  ice 
particles  consisting  of  a  sequence  of  measurements  ,  k*',-  -h/  vhere 

*  *(.*,,  XJ  4  C; 


The  measurement  error  is  represented  by  Pi  . 

In  passing  we  mention  that  other  kinds  of  measurements  of  ice  motion 
have  been  made.  Hunkins  (1967)  and  McPhee  (1978)  for  instance  inferred  the 
ice  velocity  using  current  meters  suspended  from  the  ice  into  what  they 
assumed  to  be  essentially  a  static  ocean.  These  data  were  used  to  study 
motion  on  the  time  scales  of  hours.  An  attempt  to  study  motions  on  much 
shorter  time  scales  using  accelerometers  was  made  by  Craig  (1972). 

Many  techniques  have  been  used  to  measure  ice  trajectories.  These 
are  summarized  in  Table  I.  Each  technique  has  its  good  and  bad  features. 
Over  the  past  decade  most  data  have  been  obtained  by  satellite  positioninr 
which  works  somewhat  as  follows.  Suppose  a  stable  frequency  ij  is  trans¬ 
mitted  by  a  device  on  the  ice.  The  signal  received  at  the  rapidly  moving 
satellite  will  have  a  frequency  /  which  has  been  shifted  by  the  Doppler 
effect:  -f- «  ^  *  A-f  .  The  Doppler  shift  i/  is  related  to  the  rate  of 
change  of  distance  between  the  device  on  the  ice  and  the  satellite 

jL  -  X;f»l  .  The  received  frequency  f  is  measured  at  several 

times.  If  the  satellite  coordinates  are  known  at  these  times,  then  each 
measurement  produces  an  equation  with  unknowns  and  Kct  •  Generally 
several  measurements  are  made  during  the  10  to  20  minutes  it  takes  for  the 
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satellite  to  pass  by.  The  several  equations  are  solved  simultaneously  for 
l  and  ye;u  .  Precise  positioning  requires  a  Btable  transmit  frequency, 
precise  measurement  of  the  received  frequency,  and  precise  knowledge  of  the 
satellite  coordinates.  With  care,  errors  can  be  controlled  to  the  order  of 
tens  of  meters,  ss  in  the  best  uses  of  the  navy  Navigation  Satellite  System. 

Satellite  systems  which  serve  primarily  to  relay  data  from  automatic 
data  platforms  (or  buoys)  to  data  processing  centers  also  determine  the 
location  of  the  transmitting  data  platform,  using  the  same  Doppler  positioning 
principle.  The  ARGOS  system  currently  on  the  NOAA-B  satellite  can  relay 
data  from  up  to  30  sensors  and  determine  platform  locations  to  an  accuracy 
of  a  few  hundred  meters  about  ten  times  per  day.  Fully  automatic  platforms 
cost  in  the  neighborhood  of  $6,000  with  additional  costs  depending  on  the 
desired  sensors. 

Satellite  imaging  systems  can  also  be  used  to  measure  ice  motions,  pro¬ 
vided  that  some  features  on  the  ice  can  be  identified  in  a  sequence  of  images. 
Because  of  their  all  season  and  all  weather  capability  and  good  resolution, 
imaging  radar  systems  will  probably  be  best.  The  basic  angular  resolution 
of  these  systems  is  approximately  the  ratio  of  the  wavelength  of  the  radar 
signal  to  the  diameter  of  the  antenna.  To  achieve  an  angular  resolution  of 
TO-5  (10  meters  at  a  range  of  1,000  kilometers)  with  a  wavelength  of  25  cm 
(1.2  GHz)  requires  an  antenna  25  km  in  diameter.  Although  such  large  antennas 
cannot  be  constructed  in  space,  it  is  possible  to  synthesize  large  antennas 
by  using  data  from  several  points  along  the  satellite's  orbit.  The  determi¬ 
nate  ens  of  the  geographical  position  of  an  ice  feature  with  the  SEASAT  data 
contained  errors  of  up  to  3  kilometers  (see  Hall  and  Rothrock,  198l).  Unless 
these  errors  can  be  reduced  the  data  are  not  particularly  valuable  for 
measuring  the  displacements  over  intervals  of  a  few  days  or  less.  However, 
the  errors  are  highly  correlated  in  space  and  are  essentially  eliminated  in 
estimates  of  the  spatial  variability  of  the  ice  motion. 

The  attractive  feature  of  imaging  radar  systems  is  their  potential  to 
sample  densely  in  space.  Hall  and  Rothrock' s  work  suggests  that  it  will  be 
possible  to  track  roughly  one  feature  per  square  kilometer,  which  will 
resolve  most  of  the  spatial  structure  of  the  field  of  motion.  Techniques 
for  extracting  data  from  the  images  or  from  the  raw  data  are  still  rather 
primitive.  No  doubt  satisfactory  automated  techniques  for  identifying  and 
tracking  features  will  be  developed  when  the  need  arises.  At  present  there 
is  no  imaging  radar  system  in  space.  Planning  is  underway  for  a  system  to 
be  in  operation  perhaps  by  1985 . 
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5.1.1.  Sources  of  data 

The  earliest  ice  motion  data  are  the  trajectories  of  ships  beset  in 
the  ice.  These  are  followed  by  the  trajectories  of  numerous  Soviet  and  U.S. 
drifting  research  stations,  and  more  recently  by  the  trajectories  of  auto¬ 
matic  data  buoys.  The  following  list  of  data  sources  is  by  no  means  complete. 


Experiment 

Reference 

Description 

T-3 

Hunkins  et  al.,  1971 

Monthly  positions  May  1962- 
December  1970 

Beset  ships  and 
drifting  stations 
before  1970 

Hastings,  1971 

A  chart  showing  most  prior 
trajectories,  one  point  per 
month;  includes:  Fram,  Jeanette, 
Maud,  Sadko,  Sedov,  Tecetthoff. 
British  Trans-arctic  Ext., 
Alpha,  Arlis  I,  Arlis  II, 
Charlie,  North  Pole  1-20,  T-3. 

AIDJEX  1972 

AIDJEX  Staff,  1972 

100  km  triangle,  Beaufort  Sea, 
March  and  April  1972. 

A1BJEX 

1975-76 

Thorndike  and  Cheung,  1977 

Manned  camp  positions  and 
velocities  tabulated  at  6  hr 
intervals,  daily  buoy  positions 

1976-77  buoys 

Thorndike  and  Cheung,  1977 

Tabulated  daily  positions  of 
Beaufort  Sea  buoys. 

1979  buoys 

Thorndike  and  Colony,  1980 

Tabulated  daily  positions, 

25  buoys,  analyses  of  surface 
pressure  plotted  daily. 

1980  buoys 

Thorndike  and  Colony,  1981 

As  above. 

Fram  I ,  1979 

Hunkins  et  al.,  1979 

March-May  1979;  81»°li,  9°’"’ 

LOREX  79 

Popelar  et  al.,  1981 

Three  stations,  100  km  spacing 
April,  May  1979;  88°-90°Ii. 

5.2  The  general  circulation 

The  main  features  of  the  long  term  circulation,  Figure  2,  are  the  clock- 
vise  circulation  in  the  Beaufort  Sea — and  the  motion  of  ice  from  the  Siberian 
coasts  across  the  North  Pole  and  through  the  Greenland-Spitsbergen  passage. 
Time  honored  nomenclature  for  these  features  are  the  Beaufort  Gyre  and  the 
Transpolar  Drift  Stream.  Some  handy  numbers  for  these  long  term  features  ar„* 

center  of  Beaufort  Gyre:  80°N,  155°W,  half  way  between  Pt.  Barrow, 
Alaska  and  the  North  Pole 


time  to  make  1  circuit:  5  years 

time  to  traverse  Transpolar  Drift  Stream:  3  years 

area  flux  through  Oreenland-Spitsbergen  Passage: 

300  km*/day  or  about  205*  of  the  area  of  the  basin  per  year. 

This  pattern  of  motion  exists  only  as  an  average  over  several  years. 

On  shorter  time  scales  there  are  departures  from  the  long  term  pattern.  Con¬ 
sider  for  example  the  trajectories  plotted  in  figure  3  for  the  year  1979  and 

I 

figure  for  1980.  The  trajectories  are  characteristically  meandering  and 
convoluted  showing  that  on  daily  and  monthly  time  scales  the  ice  motion 
differs  markedly  from  the  long  term  mean.  Notice  the  major  anomaly  which 
occurred  in  the  summer  of  1980  when  the  motion  of  several  buoys  for  several 
months  was  in  the  opposite  sense  from  the  long  term  mean.  These  departures 
of  the  actual  motion  from  the  long  term  pattern  in  most  cases  represent  the 
response  of  the  ice  to  the  passage  of  atmospheric  systems. 

If  we  ignore  the  forces  and  examine  only  the  ice  motion,  the  departures 
from  the  mean  circulation  appear  aperiodic  and  chaotic.  The  departures  can 
be  thought  of  as  random  but  they  are  not  without  structure.  Our  objective 
now  is  to  clarify  this  structure. 


5.3  Ice  velocity 

The  ice  velocity  field  can  be  defined  by  the  relationship 

*M>  *  l-  . 
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Here  the  particle  label  /  plays  no  role. 

This  definition  is  meaningful  only  if  the  limit  exists.  From  a  practical 
point  of  view  it  is  useful  only  if  the  limit  exists  and  is  approached  when 
the  interval  R  decreases  to  the  time  interval  T  between  observations. 

In  figure  5  several  sets  of  observations  of  ice  motion  are  plotted, 
showing  the  variation  of  one  coordinate  of  position  versus  time.  Successive 
data  sets  divide  the  sampling  interval  't  by  15  and  improve  space  resolution 
by  the  same  factor,  giving  a  sequence  of  closer  and  closer  perspectives  on 
the  motion.  By  constructing  velocity  estimates,  u^ft)  *  )]/ki 

for  k  decreasing  to  1,  one  can  examine  the  limiting  process.  In  the  first 
two  figures  of  the  sequence,  u.^  continues  to  change  appreciably  for  small  k  . 
In  the  last  two  figures,  HfcU)  at  most  times  "i  becomes  almost  independent 
of  k  for  small  k  ,  implying  that  on  these  time  scales t  X  <  15  minutes, 
the  ice  indeed  possesses  a  velocity. 
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Convergence  of  [*(*♦*,$  -  to  a  definite  limit  implies  that 

the  increment  is  proportional  to  x  for  small  x  . 

For  random  processes,  a  useful  condition  for  convergence  is  mean  squared 
convergence  in  which  the  variance  of  the  increment  becomes  proportional 
to  X*  for  small  X  .  This  can  be  readily  tested  for  ice  trajectories  by 
plotting  -*(*♦#))''  versus  t  on  a  log-log  plot,  as  is  done  in 

figure  l  .  The  process  is  differentiable  if  the  graph  has  slope  of  2  for 
small  X  . 

From  a  practical  point  of  view  the  process  ceases  to  be  differentiable 
at  measuremen  intervals  X  for  which  the  slope  of  the  graph  departs  appre¬ 
ciably  from  2.  The  evidence  in  figures  5  and  i  implies  that  the  ice  does 
have  a  velocity  and  that  it  can  be  resolved  with  a  sampling  interval  of 
about  three  hours. 

Time  averaged  velocities  can  be  defined  without  reference  to  the  limit¬ 
ing  process.  Let  u.('b,T(*/)  ^he  time  averaged  velocity  at  ■£  . 


*  tJ 

f  q-T/2 


_L  f 


This  quantity  MX*)  and  its  properties  depend  on  the  duration  T  of  the 
time  averaging.  For  example,  the  variance  of  will  in  general 

be  less  than  the  variance  of  because  the  T  average  has  suppressed 

contributions  to  the  variance  on  shorter  time  scales.  Also,  from  a  dynamical 
point  of  view,  the  equations  which  satisfies  should  involve  J  as 

a  parameter.  Different  physical  processes  may  be  responsible  for  determining 
uUiT  *)  for  different  values  of  7*. 

Typical  ice  velocities  range  from  0  to  20  cm  sec"1.  An  extreme  velocity 
of  lUO  cm  sec”1  has  been  observed.  Two  histograms  of  ice  speed  are  shown  in 
figure  7  ,  one  corresponding  to  a  full  year  of  observations,  the  other 

restricted  to  summer  observations.  The  winter  data  contained  several  periods 
of  essentially  zero  motion.  During  the  summer  the  ice  was  never  observed  to 
8top. 

In  the  following  pages  particular  importance  is  attached  to  the  tine 
and  apace  variability  of  the  ice  motion.  We  will  often  refer  to  the  variance 
of  velocity,  *(v-W  J:f\  This  quantity  has  been  evaluated  from  many 

observations j  it  varies  appreciably  with  season  and  with  location.  Table JT 
gives  estimates  of  the  velocity  variance  for  eech  of  the  buoys  shown  in  figure  3 . 
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TABLE  II 


5.3.1  Velocity  time  series 

The  ice  velocities  deduced  from  trajectories  should  be  interpreted 
from  the  Lagrangian  point  of  view  since  the  measurements  are  made  following 
the  material  particle.  If  many  trajectories  are  measured  simultaneously 
it  is  possible  to  obtain  by  Interpolation  velocity  time  series  at  fixed 
Eulerian  points.  Although  the  Lagrangian  description  is  more  directly 
measurable,  it  has  the  drawback  of  sampling  both  the  time  and  space  varia¬ 
tions  in  velocity.  The  differences  between  the  Lagrangian  and  Eulerian 
points  of  view  is  not  of  great  practical  importance  in  the  motion  of  sea 
ice.  This  is  because  most  of  the  variability  in  the  ice  motion  is  driven 
by  the  wind,  and  the  space  patterns  in  the  wind  field  move  across  the  basin 
so  much  faster  than  the  ice  moves,  that  all  points,  Lagrangian  or  Eulerian, 
experience  essentially  identical  stochastic  forcing.  Consequently,  Lagrangian 
and  Eulerian  time  series  of  sea  ice  velocities  look  Just  about  the  same. 

A  typical  velocity  time  series,  taken  from  an  ice  station  trajectory, 
say,  will  have  a  mean  velocity  usually  less  than  3  cm  sec"1,  and  it  may  show 
a  trend.  When  these  effects  are  removed  from  the  time  series  what  remains 
are  the  fluctuations  with  time  scales  longer  than  the  sampling  interval 
and  shorter  than  the  length  of  record  1~  .  The  AIDJEX  position  data  for 
instance,  serve  to  resolve  fluctuations  on  time  scales  from  a  few  hours 
up  to  a  few  months,  (  X  »  1  hour,  T *  1  year).  Over  this  range  of  time 
scales  the  ice  velocity  has  a  power  spectral  density  as  sketched  in  figure  8  . 

The  power  spectral  density  is  plotted  for  positive  and  negative  frequen¬ 
cies  corresponding  to  counter-clockwise  and  clockwise  rotations  of  the 
velocity  vector.  At  the  end  of  this  section  an  algorithm  is  given  for  cal¬ 
culating  the  spectra  of  vector  time  series.  We  interpret  the  velocity 
vector  time  series  *)  as  a  complex  time  series  with  the  Fourier 

decomposition 


The  power  spectral  density  is  the  real  function  Sew)  s  ,  defined 

in  the  frequency  range  <  40  <  "at  * 

The  integral  of  the  spectral  density  over  this  frequency  range  is  the 
variance  a*  ,  and  the  integral  over  any  frequency  band  is  the  part 
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of  the  total  variance  contributed  by  fluctuations  with  frequencies  in  that 
band.  A  useful  way  to  summarize  the  information  contained  in  the  spectrum 
is  to  state  the  fraction  of  the  total  variance  coming  from  frequencies 
greater  than  a  certain  value.  We  find; 

Table  III  Fraction  of  total  variance  from  frequencies  exceeding 
frequency:  1  cycle  per  month  2  1*  15  30  60 

period:  1  month  2  weeks  1  week  2  days  1  day  12  hours 

58#  U5J*  3l*5S  12%  1%  1% 


The  equation  of  motion  for  sea  ice  balances  the  ice  acceleration  against 
air  stress,  water  stress,  pressure  gradient  forces  due  to  the  sloping  sea 
surface,  and  internal  ice  stress  gradients.  While  it  is  not  the  purpose 
here  to  examine  the  ice  dynamics,  it  is  useful  to  relate  features  of  the 
ice  velocity  spectrum  to  these  driving  forces.  In  the  central  basin  about 
75#  of  the  variance  of  the  ice  velocity  can  be  explained  by  the  local  geo- 
strophic  wind  (Thorndike  and  Colony,  1981).  In  fact,  the  ice  velocity  fluctua¬ 
tions  are  roughly  proportional  to  the  local  wind  fluctuations.  This  implies 
that  the  spectra  of  the  ice  velocity  and  of  the  wind  should  have  approximately 
the  same  shapes  (see  figure  ?  ). 

The  water  stress  depends  on  the  difference  in  velocity  between  the  ice 
and  the  upper  ocean.  If  the  ocean  is  at  res1"-  the  water  stress  is  simply  a 
drag  opposing  the  ice  velocity  and  its  only  effect  is  to  reduce  the  ice 

response  near  the  inertial  frequency. - 

If  the  ocean  is  in  motion,  the  ice  will  be  carried  along  with  it,  and  should 
acquire  spectral  traits  similar  to  those  of  the  ocean.  Unfortunately  the 
spectral  signature  of  motion  in  the  Arctic  Ocean  is  poorly  known.  The  long 
term  circulation  of  the  upper  ocean  appears  to  be  similar  to  that  of  the  ice, 
with  a  clockwise  circulation  in  the  Beaufort  Sea  and  a  transpolar  current 
flowing  from  Siberia  through  the  Greenland-Spitsbergen  passage.  It  probably 
is  not  productive  to  ask  whether  the  ice  drives  the  long  term  ocean  circula¬ 
tion  or  vice  versa  because,  in  the  long  term,  the  ice  should  be  thought  of 
as  part  of  the  upper  ocean.  In  any  case  the  long  term  ocean  behavior  does 
not  affect  the  spectrum  we  are  considering  since  we  have  subtracted  out  the 
long  term  mean  velocity. 
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Fluctuations  in  ocean  currents  on  shorter  tine  scales  have  been  observed 
but  it  is  not  yet  possible  to  say  how  great  their  effect  on  the  ice  velocity 
nay  be.  Briefly  they  are: 

1.  Monthly  variations  are  apparent  in  the  currents  deduced  from  vind  and 
ice  notion  observations  by  Thorndike  and  Colony,  1961. 

2.  Manley's  (1981)  study  of .subsurface  eddies  during  AIDJEX  described 
nunerous  features  with  velocities  of  about  50  cm  sec-1,  a  length  scale  of 
10  kilometers,  and  a  depth  of  50-200  meters.  He  found  no  expression  of 
these  structures  in  the  ice  motion. 

3.  Inertial  oscillations.  The  balance  between  the  ice  acceleration  and 
the  Coriolis  force  leads  to  oscillatory  ice  notion  with  a  period  which  varies 
from  12.77  hours  at  70°N  to  12  hours  at  90°N.  Inertial  motions  are  always 
clockwise— hence  the  negative  value  for  the  frequency— in  the  Northern  Hemis¬ 
phere  because  the  Coriolis  force  always  accelerates  the  velocity  to  the  right. 
Inertial  motions  in  sea  ice  were  first  described  by  Hunkins  (1967)  and  have 
received  subsequent  study  by  McPhee  (1979)  and  Colony  and  Thorndike  (19So). 

The  amplitudes  can  reach  0.20  m  sec"1  during  summer  when  the  ice  pack  is 
comparatively  loose.  Their  effect  on  the  summer  ice  velocity  spectrum  is 
striking  (see  figure  0). 

1*.  Tides.  Tidal  currents  in  the  central  basin  are  small  because  of  the 
great  depth.  Theoretical  estimates  are  in  the  range  of  1-2  cm  sec”1.  Over 
the  shallow  continental  shelves  the  amplitudes  are  predicted  to  be  at  least 
an  order  of  magnitude  greater  (Kowalik  and  Untersteiner,  197f  ).  These 
theoretical  estimates  are  for  the  lunar  semi-diurnal  tide.  Evidence  from 
tide  gauges  around  the  basin  summarized  by  Sverdrup  (1926)  imply  that  the 
lunar  semi-diurnal  (period  12.1*7  hr)  and  the  solar  semi-diurnal  (12  hr)  tides 
are  the  dominant  tidal  constituents  in  the  Arctic  Ocean. 

The  prediction  of  a  large  amplitude  of  the  tidal  current  over  the  shelf 
is  confirmed  by  Sverdrup's  observations  from  the  Maud;  Nansen  may 

also  have  observed  tidal  motion  in  the  ice  pack  surrounding  the  Fram.  Because 
the  tidal  and  inertial  periods  are  so  close  it  may  be  difficult  to  diagnose 
observed  motions  correctly.  Nevertheless  there  are  several  differences 
between  the  two  kinds  of  motions  which  can  sometimes  be  used  to  distinguish 
them.  First  the  tidal  vector  traces  out  an  ellipse  during  one  period.  An 
ellipse  can  be  viewed  as  the  sum  of  a  clockwise  circle  and  counter-clockwise 
circle.  Thus,  unless  the  tidal  ellipse  should  happen  to  be  exactly  a 
clockwise  circle,  it  should  have  some  expression  on  the  counter-clockwise 
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tide  of  the  spectrum.  Inertial  motion  is  strictly  clockwise.  Second*  at 
a  fixed  point,  tidal  motions  should  have  a  fixed  phase.  Inertial  motion 
on  tne  other  hand,  acquires  a  new  phase  every  time  the  ice  receives  a 
sudden  impulse  of  momentum.  Third,  the  phase  of  tidal  motion  should  vary 
smoothly  and  slowly  in  space.  Inertial  motions  at  different  points  may 
have  no  fixed  phase  relationships. 

The  motion  described  hy  Sverdrup  can  hardly  be  mistaken  for  inertial 
motion.  The  smooth  variations  in  phase  as  the  Maud  moved  from  Vrangel  Island  to 
the  New  Siberian  Islands,  and  the  clearly  elliptical  cycle  described  by  the 
measured  velocity  vector  are  not  compatible  with  inertial  motions. 

Most  ice  trajectories  which  have  been  analyzed  in  the  western  literature 
have  been  over  the  deep  basin  where  the  tidal  motion  is  small.  Careful  obser¬ 
vations  and  analysis  would  be  required  to  detect  a  tid-.  component  of  order 
1  cm  sec”1  since  it  would  be  mixed  with  an  inertial  component  which  is  often 
much  larger  and  with  a  rich  spectrum  of  other  types  of  motions. 

The  Soviet  literature  contains  many  references  to  tidal  motion  in  sea 
ice.  Doronin  and  Khiesin  (1977)  and  Zubov  (19^3)  each  devote  several  pages 
to  the  subject. 

An  important  consequence  of  the  tidal  motion  is  the  associated  cycle 
of  convergence  and  divergence  caused  by  the  difference  in  phase  of  the  tidal 
cycle  at  different  points.  Periodic  opening  and  closing  of  the  ice  in  the 
shallow  seas  has  an  effect  on  the  heat  exchange  between  the  atmosphere  and 
the  ocean  and  on  the  rate  of  ice  production.  The  theoretical  calculations 
of  Kowalik  and  Untersteiner  indicate  maximum  divergence  rates  in  the  shallow 
seas  exceeding  10” ‘sec”1  which  is  enough  to  produce  one  percent  opening 
during  the  tidal  cycle.  Their  theoretical  estimates  of  the  divergence  rate 
associated  with  tides  over  the  deep  ocean  are  three  orders  of  magnitude  smaller. 

Calculation  of  cross  power  spectral  density  for  two  two-dimensional  vector 
time  series. 

Given:  two  discrete  complex  time  series  Vk  ,  with  sampling 

interval  r  . 

Step  1.  Select  M  ,  the  number  of  lags. 

Step  2.  Remove  the  mean  and  trend  from  each  time  series. 

Step  3.  Compute  cross  covariances 
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Define  spectral  window 
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Step  5.  Calculate  spectrum  b  and  frequency  tj 

*  -ZwJLkFT 

S,  =  T  L  rl  6 
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k=  -M, 


M 


A. 

2riz 


5.3.2  Velocity  correlations  in  tine. 

An  alternative  characterization  of  a  randor  ^unction  is  its  autocorre¬ 
lation  function.  The  discussion  in  the  previous  section  of  the  power  spectral 
density  of  the  time  function  at  a  fixed  or  moving  point  could  have  been 

given  in  terms  of  the  autocorrelation  function 


R  (?) 


E  ull+-c)u.*(Q 
£  uCt)  vStt) 


r 


The  two  functions  S(^)  and  R(z)  are  Fourier  transforms  of  each  other  and 
therefore  cojfjtain  equivalent  information.  Which  description  is  the..nore  useful 
depends  on  the  application.  The  spectrum  is  useful  for  distinguishing  physical 
processes  with  distinct  characteristic  frequencies;  for  example  separating 
the  free  inertial  oscillations  from  the  wind  forced  motion.  On  the  other 
hand  the  autocorrelation  function  is  more  useful  for  questions  related  to 
prediction  or  experiment  design.  For  example  the  question:  "How  well  can 
tomorrow’s  ice  velocity  be  predicted  on  the  basis  of  today’s  velocity,”  has 
an  answer  involving  the  autocorrelation  function  £(t)  evaluated  at  t  *  1  day. 

The  autocorrelation  function  is  complex.  Its  real  part  contains  informa¬ 
tion  about  the  lagged  correlations  of  the  pL  component  of  veloc.  y  with  itself 


and  the  ^  component  with  itself.  The  imaginary  part  has  information  about 
the  lagged  correlations  between  the  A  and  ^  velocity  components.  If  the 
velocity  of  a  piece  of  ice  undergoing  an  inertial  oscillation  has  a  positive 
A  component  at  a  certain  time,  it  will  have  a  negative  g  component  three 
hours  later.  The  imaginary  part  of  the  correlation  should  have  negative 
peaks  at  3  hours,  15  hours,  etc. 

The  autocorrelation  function  for  the  velocity  time  series  at  the  AIDJEX 
ice  station  Caribou  is  shown  in  figure  H  .  The  autocorrelation  function 
estimated  by  the  1979  buoy  data  for  the  central  basin  is  shown  in  figure  /Z  . 
The  real  part  of  the  correlation  falls  to  about  0.7  after  one  day,  0.1*  after 
two  days,  and  decreases  slowly  at  longer  lags.  The  large  correlation  at  lag 
one  day  indicates  that  persistence  (the  forecast  strategy  which  predicts 
that  the  future  will  be  the  sane  as  the  present)  will  have  sone  skill  for 
one  day  forecasts. 

As  expected  the  autocorrelation  fron  Caribou  has  negative  peaks  in  the 
imaginary  part  at  3  hours,  15  hours,  etc.  Generally  though  the  imaginary 
part  is  small.  If  the  inertial  motions tare  not  of  interest  in  a  particular 
application,  the  imaginary  part  of  the  autocorrelation  function  can  safely 
be  ignored.  This  is  equivalent  to  treating  the  two  velocity  components  as 
independent  time  series. 

5. 3. 2.1  Application  of  the  time  autocorrelation  function. 

It  is  desired  to  estimate  the  ice  velocity  u.  at  time  given  observa¬ 
tions  ^  at  times  *(,•••.»  N  .  Suppose  the  observations  have  zero 

mean  random  errors  -  uCtj)  which  have  covariance 

-  cr-^  . 

The  delta  function  expresses  the  independence  of  errors  at  different  times. 
Finally  suppose  the  errors  are  independent  of  the  actual  velocities.  Choose 
an  estimator  It  which  is  a  linear  sum  of  the  observations. 

£  -  f. 

yi  *  • 

It  is  desired  to  find  the  complex  constants  which  give  the  best  estimate 

a  in  the  sense  that  E(u-u)(u-u.)*  is  minimized.  This  is 

accomplished  by  differentiating  f  with  respect  to  the  real  and  imaginary 
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parts  of  each  and  equating  to  zero.  The  result  is  the  system  of 
complex  equations  for  the  unknowns 

■  (fv*/ 

A  further  reduction  is  achieved  by  noting  that 

f.j.;  .  ev*  -  «■*£,•«  E  *iu‘  ■  *V*  • 

Use  of  the  matrix  notation  <?  '•  ,  A-  t#ji  ,  ^  j 

gives  the  compact  expressions 

£  =  -at2 

M#A  m  P  . 

Furthermore  the  estimation  error  is 

F  «  $  -  . 


w 


2 


With  these  two  expressions  we  can  answer  a  number  of  practical  problems. 
1.  What  sampling  rate  is  required  to  insure  good  two  point  interpolation 
in  time?  For  given  sampling  interval  r  and  measurement  error  variance  <rx  , 
we  write  eq.  /  as 


The  solutions  for  the  o<^  ,  which  in  this  case  are  complex  conjugates  for 

reasons  of  symmetry,  and  for  F  are  sketched  in  figure  /5 .  The  solution 
technique  extends  trivially  to  interpolation  involving  more  points. 

The  optimal  choice  of  the  weights  and  ***  is  not  «*,  *  **2  ■  &  .  This 

linear  interpolation  between  the  two  data  points  is  the  most  natural  scheme, 
and  it  is  quite  good  for  small  T.  ,  and  small  cr  ,  but  it  is  not  the 
optimal  interpolation.  Especially  at  large  r  ,  a  smaller  error  variance 
is  achieved  by  giving  less  weight  to  the  observations.  For  very  large  t 
the  optimal  estimate  is  simply  .  (.•‘•'•V0,)  since,  in  this  case  the 
two  observations  are  so  removed  in  time  as  to  have  no  correlation  with 
the  desired  velocity. 
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The  measurement  error  tr  affects  and  F  slightly.  The  dependence 
of  F  on  the  sampling  interval  X  is  quite  strong.  It  appears  that  if 
one  needs  to  know  each  velocity  component  to  i  3  cm  sec”'  ,  F  ■ £  *1  *  = 

18  cm  sec”',  a  sampling  interval  of  about  six  hours  is  required  even  with 
perfect  measurements. 


2.  What  errors  are  expected  in  forecasts  based  only  on  the  present  velocity? 


Here  we  obtain 


/  +  tr’/f 1 


"/  y  o- 1  /f * 


The  forecast  skill  is  plotted  for  several  values  of  initial  error  in  figure  14 
The  prediction  error  grows  with  time,  approaching  the  velocity  variance 
for  large  7  .  From  the  sketch  we  see  that  an  optimal  3  hour  forecast 
captures  about  75/S  of  the  velocity  variance,  an  optimal  one  day  forecast 
captures  about  35^  of  the  variance. 

The  velocity  was  assumed  to  have  zero  mean  in  this  analysis  and  in 
constructing  the  autocorrelation  function  (figures H/i)  the  long  tern  mean 
was  first  removed  from  the  data.  In  an  actual  application  to  a  region  where 
the  mean  is  known  the  best  forecast  would  be 

U  -  u.  +  (*  -  * )  • 


Note  that  the  best  forecast  is  not  simple  persistence  (  ®<  *  1).  For  long 

( 

forecasts,  the  best  estimate  is  the  meanA  ,as  equation  o  correctly  indicates 
5. 3*2.1.  Acceleration 

The  time  autocorrelation  function  for  velocity,  figures  ll  and  12  ,  behave 
like  /-**/£  for  small  'l  .  This  implies  that  the  velocity  is  not  tine 

differentiable  in  the  mean  square  sense.  Attempts  to  measure  accelerations 
precisely  by  sampling  over  shorter  and  shorter  time  intervals  may  merely 
result  in  larger  and  larger  estimates  of  acceleration.  Of  course,  time 
averaged  accelerations  exist  and  have  finite  variance. 

g  ‘  .  v.r  {4  M*  )  •'  ( ' '  '  * ' 17 $ 

which  for  small  is  2f*"/cT  vith 
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It  is  frequently  claimed  that  the  inertial  forces  due  to  acceleration 
of  the  ice  are  small  compared  to  other  forces  acting  on  the  ice.  To  support 
this  claim  a  typical  average  acceleration  05.  is  compared  for  instance,  to  a 
typical  value  of  the  Coriolis  acceleration  ‘fo*  *  •  |0  ***••*■' 

Using  the  above  estimate  for  (£_  ,  and  taking  ^  »  teo  e/«»*  ,  we  find 

<  £0^  provided  T  >  2 'ID5  .  Thus,  roughly  a  half  hour  time 

average  is  required  to  bring  the  acceleration  down  to  the  level  of  other 
terms  in  tne  momentum  balance. 

Velocity  correlations  in  snace 


In  earlier  sections  the  emphasis  has  been  on  the  structure  of  the  complex 
function  U.  of  the  real  argument  4  ,  at  a  fixed  or  moving  point.  In  this 
section,  the  emphasis  is  on  U.  as  a  function  of  position  for  fixed  time. 
Certain  results  are  more  easily  expressed  if  U.  and  ^  are  regarded  as  two 
dimensional  vectors  rather  than  as  complex  numbers.  The  autocorrelation 
function 


£ul(tOv|*;  E  v(*,)W 9t)J 


is  a  2  x  2  matrix  involving  possibly  four  different  functions,  whereas  the 
time  autocorrelation  function  involves  Just  two. 

The  structure  of  the  spatial  autocorrelation  matrix  for  sea  ice  veloci¬ 
ties  reveals  some  of  the  properties  of  the  velocity  field.  To  establish 
the  basic  ideas  we  fir«t  determine  the  structure  of  the  autocorrelation  matrix 
for  a  homogeneous,  isotropic,  two  dimensional  velocity  field,  borrowing  from 
Kolmogroff's  (19^1 )  classic  discussion  of  three  dimensional  turbulence. 

A  random  field  is  said  to  be  homogeneous  if  its  statistical  properties 
are  invariant  under  a  translation  of  coordinates,  and  isotroric  if  they  are 
invariant  under  rotations  and  reflections  of  coordinates.  In  particular  the 
autocorrelation  matrix  for  a  homegeneous  isotropic  field  must  satisfy 


ft-  )  for  any  ft  ,  and 

for  any  orthonormal  matrix  ^1  .  Multiplication  by  an  orthonormal  matrix 


accomplishes  an  arbitrary  rotation  or  reflection  of  coordinates.  Note 
that  the  change  of  coordinates  is  applied  to  both  the  position  and  the 
velocity  vectors.  These  definitions  make  precise  the  notions  that  the 
statistics  do  not  depend  on  position  or  direction.  By  choosing  the  trans¬ 
lation  vector  a  or  the  rotation  and  reflection  matrix  M  suitably,  ve  can 
see  how  these  definitions  of  homogeneity  and  isotropy  constrain  the  form  of 
the  autocorrelation  matrix. 

First,  *i)  depends  only  on  the  vector  difference  .  To 

see  this  take  a  •  -*<  .  Note  that  this  does  not  imply  that  depends  only 

on  the  distance 

Second,  k  is  symmetric.  This  requires  using  «.•"£'  »  ^  '  and 

o.  t  %  *  to  obtain  this  sequence  of  equalities: 

0m  m 

Third,  /V/o,  0/V(*,  •  H6re  use  the  reflection  invariance  in  the 

definition  of  isotropy  by  taking 


M 


■  (:  -.) 


which  changes  the  sign  of  the  £  axis.  Then 

furo.  >)*(?,*)  *  f 

which  is  only  true  when  £u(o,t>)  rO.  / 

Consequently  the  autocorrelation  function  evaluated  at  7,vl °<e)  •  7t 


has  the  form 


E  uK*}**tT»op 


involving  only  two  function  of  distance.  Byl*)  is  the  correlation  between 
the  components  of  velocity  parallel  to  the  line  Joining  two  points  separated 
hy  a  distance  r.  is  the  correlation  between  the  components  of 


velocity  perpendicular  to  that  line. 

More  generally  if  *  and  X,  are  arbitrary  points  with 

OM&HVJO/r  and  »  then 


r**  (»,-*»/*  ’If ‘A 
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which  is  the  general  form  for  the  autocorrelation  matrix  for  a  homogeneous , 
isotropic  two  dimensional  velocity  field.  It  involves  two  functions  Q\\ 
and  %  which  we  will  estimate  below  from  observations.  Further  it  involves 
the  separation  *♦*  and  the  orientation  B  of  the  two  points  *■  and  at 
which  the  autocorrelation  is  being  evaluated.  It  may  seem  strange  that  £ 
enters  the  statistics  for  an  isotropic  field.  The  condition  of  isotropy 
does  not  imply  that  fu(  6,o )U.L<,e)  equals  £u.tt,o)  u(t>,-r)  however.  A 

rotation  of  coordinates  affects  both  the  velocity  and  the  position  vectors. 

A  valid  statement  is  £ u(c,o)u.(r,o)  *  £  v(o,  o)  v(c,  r)  See  figure  15. 

There  is  another  constraint  on  the  spatial  correlation  function.  Let 
\xL  3  (u-MvJ  be  the  velocities  at  any  set  of  points  %i  ,  >*),  and  let 
jf\  ^  be  arbitrary  real  numbers.  Then  the  linear  sum 

S'  %  u«.  7, 

»*/ 

has  a  variance  which  can  be  expressed  in  terms  of  \  .  The  added  constraint 
is  that  yar  ($}^o  for  all  choices  of  •  This  rtquires  that 

the  matrix  consisting  of  the  correlationSof  the  velocity  component  at  any 
set  of  points  [xjraust  be  positive  definite  (all  positive  eigenvalues),  and 
is  equivalent  to  the  condition  that  the  two  dimensional  spatial  power  spectrum 
of  velocity  be  positive  for  all  pairs  of  wave  numbers  (**>*},) 

Contradictions  can  arise  if  correlation  functions  are  used  which  do  not 
satisfy  this  constraint.  For  example,  suppose  we  wish  to  consider  a  velocity 
field  for  which  Bill®)’  I,  G>j(2L)  *  °*5  .  The 

attempt  to  evaluate  the  variance  of  8*  -  /.y  *  o.yu(U,o)  in  terms 

of  these  correlations  produces  a  negative  result.  Since  variance  is  intrinsi¬ 
cally  positive,  we  must  conclude  that  no  velocity  field  could  have  the  presumed 
correlations.. 
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As  an  example  of  an  analytical  fora  for  which  satisfies  the  positive 
definite  constraint  and  bears  some  resemblance  to  reality,  define 

,  GjH  -•  (  >  '  e‘ ( ,  A  -loo**  ■ 

This  example  was  constructed  by  taking  a  stream  function  5  with  the  posi- 

_  -tVa*' 

tive  definite  correlation  function  K(r)«  e  •  The  corresponding  power 

spectrum  is  the  Hankel  transform  of  *” T/A  .  This  is  4  )  , 

which  is  always  positive.  Now  define  velocity  components  u.  *  , 

v  s  "Sl/l**  Seduce  B#[r)  =•  -  and 


C'r),  - Hf_r'  . 

”c  r 1 


5. 3. 3.1  Estimates  of  the  autocorrelation  functions. 

In  order  to  estimate  the  autocorrelation  matrix,  simultaneous  velocity 
observations  are  required  at  pairs  of  points  separated  by  distances  rar.rinr 
from  roughly  1  to  10S  kilometers.  We  present  data  from  the  19"9  buoy  pro¬ 
gram  and  from  AIDJEX  which  define  the  functions  and  13^  only  for 
distances  exceeding  100  kilometers.  The  1979  data  were  first  interpolated 
to  give  velocity  estimates  at  a  uniform  grid  in  space  and  time. 


.  a*  >-  - 

The  mean  velocity  over  all  points  i.j.fc 
the  lagged  correlations  were  found  using 


,  l*  ,  js  ,  k'/j 

was  removed.  Then  for  lags  •»» 


% 


I 


N, -J  N,-» 

i  i  r 


For  these  calculations  A/,  *  7,  ^  *  <,  and  “  200.  The  results  are  sketched 
in  figure  H  for  the  lags  JL  *  1,  ...,  7,  wf  «  *»  ■  0,  and  ^  *  n  *  0, 

7n  «  1,...,  1*.  The  results  show  that  F  uic,o)u(^,  c)  and  Evit,r'  »(f,r) 

are  similar  functions  of  r  ,  and  that  F  v(  e,t)  u.(r,c.)  ,  IT  u.(c  e'vfc.r )  / 

F  u(olo)v(r,oj  >  a*J  E  V(o, a;  w  ( e, r)  are  all  small.  These  observations 

are  roughly  consistent  with  the  assumption  that  the  field  is  homogeneous 
and  isotropic.  We  will  proceed  then  on  the  assumption,  supported  by  these 
Hta,  that  the  ice  has  a  homogeneous,  isotropic  velocity  field  with  the 
functions  Gj  and  (3j_  empirically  determined  from  figure  lL  . 

Correlations  for  the  AIDJEX  data  were  constructed  by  choosing  two  of 
the  ice  camps,  resolving  their  velocities  into  components  parallel  and  peroen 
dicular  to  the  line  Joining  the  camps,  and  correlating  these  components. 
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These  correlations  are  also  plotted  in  figure  ll 
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5. 3. 3. 2  Applications  of  the  space  correlation* 

The  spatial  autocorrelation  function  is  the  basis  for  sound  experiment 
design  and  data  processing.  For  instance ,  suppose  observations  of  velocity  u i, 
are  made  at  N  points  #;  with  independent  zero  mean  errors  having  vari¬ 
ance  <r*  .  Estimates  of  velocity  are  to  be  made  at  arbitrary  points  x  using 
these  observations.  What  is  the  maximum  estimation  error?  We  use  the  tools 
developed  earlier.  To  estimate  the  u.  component  of  velocity  at  *  ,  we  use 


observations 

*•*  *!iw 

and  construct 

/  Ruu(  -mt* 

Ru.v, 

R  U,H.t  •  • 

fiv.u, 

fcv.v.  +  F,  ' 

Rv,*t 

Ru.u, 

* 

'  u, 

•  • 

R  u. 


w  W 


wh^re  *  £  u.(*; )  v(  *  p  etc.  as  given  by  Eq.  f4  .  The  esti¬ 
mation  error  variance  is  given  by  Eq.  2  ,  Figure  //  shows  how  the  maximum 

estimation  error  depends  on  the  separation  between  measurement  stations  and 
on  the  measurement  accuracy.  Tb  iurve  labelled  C/V  refers  to  the  idealized 
condition  in  which  the  measurement  errors  are  zero.  Thus  it  represents  the 
interpolation  error  which  is  due  to  the  intrinsic  spatial  variability  of  the 
velocity  field. 


The  standard  deviation  of  the  ice  velocity  itself  is  about  10  cr.  sec_/ . 
A  reasonable  observational  goal  is  to  keep  the  interpolation  errors  below  2  or 
3  cm  secW .  This  can  be  achieved  if  the  raw  measurements  of  velocity  are 
good  to  about  3  cm  sec”*  and  the  grid  spacing  is  about  1*00  km.  At  larger 
grid  spacings  the  interpolation  accuracy  deteriorates  rapidly. 

5.1*  Deformation 

The  differences  in  velocity  from  place  to  place  are  responsible  for  the 
characteristic  morphology  of  the  ice.  When  the  velocity  difference  between 
neighboring  pieces  of  ice  is  such  that  they  tend  to  move  apart,  a  lead  forms 
and  widens  exposing  the  ocean  surface  to  the  atmosphere.  During  the  winter, 
ice  growth  is  rapid  over  open  leads.  If  the  motion  changes — perhaps  because 
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of  a  change  in  the  winds— to  that  these  pieces  of  ice  move  toward  each  ot  tr, 
the  open  lead  will  close,  and  any  new  ice  which  formed  there  will  need  to 
rearrange  itself  to  occupy  a  smaller  area.  Typical  mechanisms  for  this 
rearrangement  are  rafting  where  one  part  of  the  new  ice  sheet  overrides  the 
other,  and  ridging  where  the  ice  is  crushed  into  pieces  which  pile  into 
ridges  rising  a  meter  or  two  above  and  sometimes  many  meters  below  the  sur¬ 
rounding  ice.  Bidging  and  rafting  are  not  restricted  to  thin  ice.  If  the 
closing  motion  continues,  the  original  pieces  of  ice  come  in  contact  and 
one  may  override  the  other  or,  by  grinding  together,  pieces  may  break  off 
and  pile  up  and  down  to  form  a  ridge.  The  essential  ideas  here  are  that 
the  ice  accommodates  divergent  motion  by  increasing  the  area  of  open  water 
rather  than  by  stretching  and  thinning.  It  accommodates  convergent  motion 
by  reducing  the  area  of  open  water  and  by  ridging  and  rafting.  These  processes 
link  the  morphology  of  the  ice,  characterized  by  such  features  as  floes, 
ridges,  and  leads,  to  the  spatial  differences  in  the  ice  velocity.  By 
studying  the  spatial  variations  in  velocity  we  may  be  able  to  understand 
better  why  the  ice  pack  has  the  form  it  decs. 

There  is  a  second  reason  for  studying  the  deformation  of  sea  ice.  As 
the  ice  pack  deforms,  stresses  develop  within  it  which  tend  to  oppose  the 
deformation.  These  stresses  figure  into  the  local  balance  of  forces  and 
therefore  affect  the  motion  of  the  ice.  To  be  more  precise,  the  balance 
of  forces  equation  contains  terms  of  the  form  where 

is  the  stress  tensor.  The  stresses  are  in  turn  related  to  the  ice  motion  by 
a  constitutive  law 


stress  *  p( deformation ) . 

One  of  the  motivations  for  AIDJEX  was  to  investigate  the  function  F  from 
a  theoretical  point  of  view  and  by  using  observations  of  the  deformation  and 
indirect  estimates  of  the  stress  (see  Bothrock  et  al.,  1980). 

The  concepts  which  have  been  applied  to  the  study  of  the  spatial  varia¬ 
tions  in  ice  velocity  are  those  appropriate  to  the  analysis  of  the  deformation 
of  a  continuum.  The  fundamental  concepts  are  the  partial  derivatives  of 

and  with  respect  to  f  and  £  .  The  line  of  thought  is 

that  a  description  of  the  large  scale  deformation  in  terms  of  large  scale 
average  derivatives  should  give  some  idea  of  how  much  opening  and  closing 
is  going  on  on  smaller  scales  and  some  ides  of  the  state  of  stress.  For 
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example  under  a  purely  divergent  motion  only  leads  are  formed.  During  a 
purely  convergent  motion  no  leads  are  formed.  On  a  small  scale,  shearing 
along  a  crack  is  of  no  morphological  consequence.  But  large  scale  shearing 
is  often  expressed  locally  by  opening  at  some  cracks  and  closing  at  others. 

Mathematical  formulation  of  these  ideas  requires  two  steps.  Required 
first  is  an  appropriate  quantitative  description  of  the  large  scale  defor¬ 
mation  in  terns  of  observable  parameters.  The  second  requirements  is  for 
a  functional  relationship  giving  the  opening,  closing,  or  state  of  stress 
in  terms  of  the  observable  parameters. 

Current  practice  is  to  describe  the  motion  in  terms  of  the  large  scale 
strain  rate  invariants,  which  themselves  are  constructed  from  the  large  scale 
velocity  derivatives.  The  first  invariant 

r  .  ^ 

n 


quantifies  the  large  scale  divergence 


r  .-ay)1  +  )2' 

L  Ux  vr^ar 


and  the  second  invariant 
'/2 


cjuairt/fit*  +£*  rate  inj 


A  alkr^aVivT-  p<iV  »J-  i»v*r u 


on  ,e),  *»«■« 

If  I*  [Ei^E/J* 


expresses  the  total  rate  of  deformation  and 

0  •  o  <  6  <  r 

indicates  whether  the  motion  is  predominantly  divergence  ,  shear  &  ~  r/2. 

or  convergence  &  ~ W  • 

In  the  theory  presented  by  Thorndike  et  al.,  1975,  opening  and  closing 
are  assumed  to  be  known  functions  of  9  and  proportional  to  }£)  . 

A\  •  total  opening  ■  1*1  *.(*>  . 

C  *  total  closing  *  |E)  <*r  (*)  • 

The  functions  •<,  and  ^■r  which  give  the  total  opening  and  closing  in 
terms  of  the  strain  rate  invariants  may  have  a  form  somewhat  as  sketched  in 
figure  lB  .  There  is  not  much  hard  evidence  to  base  these  sketches  on.  In 
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fact  the  functional  relationships  lor  A  and  C  in  terms  of  |F|  and  6 
can  only  he  regarded  as  approximate  since  they  do  not  take  into  account  the 
actual  geometry  of  the  floes  in  the  region  being  considered.  Different 
regions  with  different  geometries  vill  respond  differently  to  similar  large 
scale  strain  rates. 

In  practical  applications  of  these  ideas  the  procedure  has  been  to 
take  measured  velocities  at  several  points  in  a  region  of  interest— perhaps 
the  three  corners  of  a  100  km  triangle.  Velocity  derivatives  are  estimated 
from  the  measured  velocities.  The  estimated  velocity  derivatives  are  com- 
bined  to  form  the  strain  rate  invariants  |F|  and  0  ,  and  these  are  used 
to  find  the  rates  of  opening  and  closing  A  and  C  ,  and  the  state  of  stress  T  . 

In  the  following  sections  we  reexamine  some  of  these  ideas,  particularly 
the  notion  of  velocity  derivatives,  the  influence  of  the  ice  pack  geometry 
on  opening  and  closing,  and  the  difficulties  in  basing  inferences  about 
these  quantities  on  small  numbers  of  measurements. 

5.^.1  Velocity  derivatives 

Since  the  ice  pack  is  made  up  of  discrete  pieces  moving  relative  to  each 
other,  the  variations  in  velocity  have  two  types: 

1)  for  any  X,  and  Xt  on  the  same  rigid  piece,  the  linear  relation 

(/>*»  t*  \  must  hold,  where  o> 

I  '  " 

o  Oc  u  ) 

is  a  constant  related  to  the  rate  of  rotation  of  the  rigid  piece. 

2)  for  X,  and  <*t  on  separate  pieces. 


u  (*,  )  -  u.(*, ) 


is  arbitrary. 


If  we  select  a  random  point  on  an  ice  floe,  the  partial  derivative 

of  velocity, 

2a  --  te. 

~*7-  l 


is  well  defined  because  for  small  enough  £  ,  and  almost 

always  lie  on  the  same  floe  and  the  linearity  of  the  first  type  of  motion 
mentioned  above  implies  the  existence  of  the  limit.  Thus  the  partial  deriva< 
tive  exists,  but  it  describes  only  the  rigid  body  rotation  of  a  floe.  It 
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■ays  nothing  about  the  relative  motion  between  floes  which  is  after  all  the 
motion  causing  the  opening  and  closing  and  the  exchange  of  momentum  between 
floes. 

To  include  the  effect  of  velocity  differences  between  floes  we  cannot 
let  £  become  arbitrarily  small.  Consider  the  velocity  increment 

where  (*,j)  is  fixed  and  R  begins  at  100  km 
and  decreases  to  the  smallest  value  R,  such  that  (*+1,, ))  and  (*.$;  lie 
on  separate  floes.  This  increment  captures  the  variation  in  velocity  from 
floe  to  floe  and  has  nothing  to  say  about  the  rigid  rotation  of  the  floe 
lies  on.  The  velocity  increment  will  generally  decrease  as  & 
approaches  ,  and  we  can  examine  the  characteristic  rate  of  decrease  by 
evaluating  the  variance  C  (  -  u.(*. j.-tOj  .  if  this  quantity 

is  proportional  to  ft4  for  small  R  ,  the  velocity  at  least  has  a  derivative 
in  the  mean  square  sense. 

The  best  data  available  for  examining  the  variance  of  velocity  increments 
for  a  range  of  spacings  were  extracted  from  SEASAT  synthetic  aperture  imaging 
radar  by  Hall  and  Rothrock.  By  comparing  two  images  separated  by  three  days, 
they  were  able  to  measure  the  displacement  of  enough  ice  features  to  resolve 
the  full  spatial  detail  In  one  dimension  of  the  ice  motion.  After  interpo¬ 
lating  to  evenly  spaced  points  #,•  their  data,  have  the  form 

,  v(*J  ]  i*l,  ,  A*  *  2 

Given  1,  . ..,  the  ^-increments  were  defined  as 

IjW*  v(*iti)  • 

The  variance  of  the  increments  was  estimated  from  the  { )}  <  * ), 
with  the  results  plotted  in  figure  1} .  The  linearity  of  the  log-log  plot 
for  small  £  supports  a  power  law  relationship  with  war  (Hi))  *  c  PT  with 
•t  in  'HU  utljUW  li«i  /•/  ^  /.J  .  Tt  appears 

then  that  the  increments  do  not  decrease  as  fast  as  ,  and  therefore 

the  velocity  does  not  have  a  derivative  in  the  mean  square  sense. 

The  statement  Just  made  for  the  variance  of  increments  can  be  recast  in 
terms  of  the  autocorrelation  function.  It  implies  that  for  small  r 

Bu(+)  .  ~  /-■£*•* 
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thus  helping  to  define  the  shape  of  0*  and  fix  at  r  *&  .  At  <r  «o  # 
and  ©x  have  the  value  unity,  are  continuous,  have  a  continuous  first 
derivative  of  aero,  hut  fail  to  have  a  second  derivative. 


5.k,g  Velocity  differences 

The  definition  of  the  derivative  as  a  limit  could  be  discarded  in  favor 
of  a  velocity  difference  over  a  finite  distance. 


and  similar  definitions  for  • 

There  is  no  question  about  the  existence  of  these  differences.  Defi¬ 
nitions  for  strain  rate  invariants  follow  naturally  by  associating  u-x.l 
with  etc.  in  the  earlier  definitions.  When  L  is  chosen  to  be  of 

order  100  ko,  these  definitions  can  be  used  to  describe  the  large  scale 
deformation. 

There  are  several  problems  with  this  approach.  First,  the  L -differences 
carry  no  information  about  smaller  scales.  Second,  no  basis  has  been  given 
for  choosing  a  particular  value  for  £.  .  Third,  measurements  of  velocity  are 
rarely  available  at  uniformly  .spaced  points.  Still  most  analyses  of  sea  ice 
strain  have  been  done  with  a  length  scale  implicitly  fixed  by  the  scale  of 
the  observations.  To  interpret  these  analyses  we  must  determine  the  role 
played  by  L  . 

Consider  the  covariances  between  the  various  ]_ -differences.  These 
follow  from  the  covariance  structure  of  the  velocity  itself.  For  example 

— 0n  l^l)  -  2  B„(r)  ♦  £?„  (r-L)J  .  (5.) 


To  interpret  this  equation,  suppose  the  [  -difference  Uj,^  is  measured 
at  two  points  separated  by  a  distance  Y  along  the  p(.  axis.  The  right  hand 
side  of  equation  S  approximates  the  second  derivative  of  at  t  . 

For  small  r  ,  G|  is  concave  down  so  the  two  L  -differences  are  positively 
correlated.  At  some  value  of  T  ,  Bg  (?)  has  an  inflection  point.  Over 
such  a  distance  the  L  -differences  become  uncorrelated;  at  longer  distances 
the  b  -differences  are  negatively  correlated.  This  example  has  been  worked 
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for  the  U*,L  component  at  points  separated  along  the  ^t-axis,  but  the 
same  procedure  can  be  used  to  correlate  any  tvo  first  velocity  differences 
at  any  tvo  points. 

It  is  apparent  from  this  and  figure  16  that  quantities  related 

to  the  -differences  are  highly  correlated  in  space  only  over  distances 
of  less  than  about  1*00  km.  Furthermore  this  statement  is  not  sensitive  to 
the  choice  of  ,  since  the  expression  S  will  change  sign  at  approximately 
the  same  T  ,  independent  of  L  .  The  correlation  length  scale  for  velocity 
by  comparison  is  roughly  twice  as  large  (see  figure  1 6). 

The  ice  velocities  are  veil  correlated  at  distances  of  up  to  100  km, 
as  evidenced  by  the  patterns  in  the  long  term  ice  motion  (figures  2-U)  and 
in  shorter  term  motions  (see  Thorndike  and  Colony,  1980  for  instance).  It 
may  seen  surpising  not  to  find  similar  patterns  in  the  strain  fields.  The 
reason  why  such  patterns  are  not  found  is  that  the  strain  field  has,  by 
equation  5,  and  figure  16,  a  correlation  length  scale  of  only  a  few  hundred 
kilometers.  Measurement  arrays  like  the  1979  and  1980  buoy  arrays  with 
a  buoy  spacing  of  roughly  500  kilometers  are  too  coarse  to  resolve  patterns 
in  the  strain  field. 

For  the  example  Just  worked  ve  can  find  the  variance  of  the  t  -difference 
by  setting  r  ‘‘  o  . 


The  symmetric  matrix  of  covariances between  pairs  of  L  -differences  at  the  same 
point  is 

( _ liiit _ V^L _ _ _ Kjf.l _ V*,  L. _ 

This  matrix  can  be  evaluated  at  L  *  ^00  km  for  example  by  reading  points 
from  the  plotted  functions  (?|  and  @L  in  figure  1^. 
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The  covariance  matrix  for  the  l  -differences  illustrates  several  properties 
of  the  vay  the  ice  deforms.  Notice  that  the  variances  of  the  parallel  dif¬ 
ferences  and  Vy<L  are  smaller  than  the  perpendicular  differences  and 

.  The  parallel  differences  relate  to  the  linear  stretching  and 
contracting  of  the  ice  pack.  The  perpendicular  differences  relate  to  shear 
and  rotation.  Apparently  the  second  type  of  motion  is  greater,  but  recall 
the  word  of  caution  given  earlier  that  the  interpretation  of  this  kinematic 
evidence  is  ambiguous.  We  cannot  say  whether  the  observed  motion  reflects 
some  property  of  the  ice  pack  or  some  characteristic  of  the  driving  forces. 

With  the  l  -differences  one  could  define  l  -strains  in  a  natural  way. 

The  statistics  in  Table  %  of  the  -strains  follow  from  the  above  matrix. 

These  values  predict,  on  the  basis  of  the  spatial  correlations  above,  that 
the  large  scale  divergence  will  typically  be  smaller  than  the  vorticity  and 
shear,  a  prediction  verified  by  many  different,  sets  of  observations. 
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Deformation  measurements 
Several  studies  have  been  made  over  the  past  decade  or  so  using  obser¬ 
vation  of  velocity  at  an  array  of  stations*  with  the  objective  of  monitoring  the 
large  scale  deformation  of  the  region  covered  by  the  array.  While  it  is  never 
stated  Just  this  way,  the  implicit  definition  of  large  scale  deformation  appears 
to  involve  spatially  averaged  derivatives  over  the  region. 


R 


where  A  is  the  area  of  the  region  sampled  by  the  array.  Similar  expressions 
define  the  other  large  scale  average  velocity  gradients.  Application  of  the 
Green-Gaus  theorem  implies 

¥  -if-  «•; »  to 

c 

where  is  the  outward  normal  to  the  perimeter  C  of  R  ,  and  i  a  unit  vector 
in  the  pi  direction. 

Note  the  difference  between  the  L  -difference  and  the  large  scale 
average  derivative.  The  former  can  be  measured  exactly  as  ; 

the  latter  requires  measurements  at  every  point  around  a  closed  curve.  In  the 
next  section  we  will  discuss  how  many  measurements  are  required  to  achieve  a 
desired  accuracy  in  the  large  scale  deformations.  But  first  we  review  some  of 
the  results  of  the  studies  Just  alluded  to. 

In  these  studies  the  deformation  estimates  were  made  by  finding  the 
linearly  varying  velocity  field  which  most  closely  fit  the  observed  velocities 
at  an  array  of  points  within  ft. 

u(*)=  u-,  +  M  * 

where  the  matrix  M  contains  the  four  large  scale  average  velocity  partial  deriva¬ 
tives.  M  and  U,  are  determined  using  a  least  squares  fit  of  the  observed 
velocities  u  and  positions  K,  .  The  strain  rate  invariants  were  then  computed 

A* 

from  the  elements  in  tf  . 

The  values,  given  in  Table  X  ,  confirm  some  of  the  results  deduced 
above  from  the  observed  spatial  correlation  functions.  For  example  it  is  clear 
that  the  vorticity  is  generally  several  times  as  large  as  the  divergence.  The 
standard  deviation  of  the  divergence  and  the  shear  are  roughly  the  same  sire, 
but  the  mean  shear  is  much  larger  than  the  mean  divergence  in  most  cases. 
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Sources? :  HiMer,  et  al.,  197U;  Colony  and  Thorndike,  1981;  Maykut,  1981. 
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The  distribution  of  the  invariant  0  is  plotted  in  figure  20  .  Note 
that  nearly  always  .  In  this  range,  the  two  principal  values 

of  the  strain  rate  tensor  have  opposite  signs.  The  figure  indicates  that 
it  is  almost  always  the  case  that  when  the  ice  pack  is  extending  in  one 
direction,  it  is  contracting  in  the  perpendicular  direction.  Ribler  et  al., 

19?1*  also  comments  on  this  for  strain  measurements  on  a  20  kilometer  scale. 

It  is  also  evident  from  Table  V  that  the  summer  values  are  usually 
somewhat  greater  than  the  winter-spring  values.  Presumably  the  ice  pack 
is  weaker  and  offers  less  resistance  to  deformation  in  the  summer.  (The 
alternative  hypotehsis,  that  the  driving  forces  are  larger  in  the  summer, 
is  not  true. ) 

The  mean  quantities  in  Table  V  produce  large  strains  over  the  course 
of  a  year.  The  year  long  deformations  at  the  AIDJEX  arrays  are  shown  in 
figure  2)  .  The  region  experienced  a  net  clockwise  rotation  of  about  35°. 

The  principal  deformation  involved  a  stretching  of  about  90%  in  the  east-west 
direction  and  a  contraction  of  about  1*0/5  in  the  north-south  direction.  The 
net  divergence  was  not  significantly  different  from  zero.  The  two  nested 
arrays  expressed  similar  deformations.  . 

Deformation  estimates  at  a  number  of  points  in  the  central  basin  for 
the  years  1979  and  1980  are  presented  in  figure  22..  The  strain  ellipses 
typically  show  large  shear  and  small  divergence.  The  only  pattern  evident 
in  the  figure  is  the  similar  alignment  of  the  major  principal  axis  of  shear 
for  the  five  points  closest  to  the  pole.  Ve  should  not  expect  to 

see  any  patterns  in  the  deformation  displayed  on  this  scale.  Recall  that 
the  spatial  correlation  function  for  /.-differences  has  a  length  scale  of 
only  about  1*00  kilometers.  The  deformations  at  points  separated  by  greater 
distances  should  evolve  more  or  less  independently. 

5.1*.l*  Interpretation  of  deformation  measurements 

A  number  of  authors  have  confronted  the  difficulties  of  describing 
the  deformation  of  this  decidedly  dis continuum.  Nye  (1973)  for  instance, 
puts  forward  a  definition  of  "strain  on  a  length  scale  L  By  first  smoothing 
the  velocity  field  using  a  kernel  of  length  scale  t*  •  he  obtains  a  new 
velocity  field  which  is  differentiable  and  for  which  the  UBual  notions  of 
deformation  based  on  partial  derivatives  are  valid.  Papers  by  Hibler  et  al., 
1971*  and  Thorndike  and  Colony,  1977  take  a  similar  point  of  view  by  attempting 
to  partition  the  ice  deformation  into  two  parts.  The  first  part  is  assumed 
to  be  associated  with  the  spatial  variation  of  the  atmospheric  and  oceanic 
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forcing  fields  and  with  the  geometry  of  the  basin,  to  have  a  characteristic 
3 

length  of  10  km,  and  to  be  differentiable.  The  second  part  is  associated 
with  the  irregular,  discrete  geometry  of  the  ice  pack.  It  has  length  scales 
less  than  10  km  and  is  characterized  by  discontinuous  variations  in  velocity. 

In  these  papers,  the  first  kind  of  variation  is  referred  to  as  the  large 
scale, underlying,  continuum,  or  linear  part  of  the  velocity  field,  while 
the  second  is  referred  to  as  the  small  scale,  local,  perturbation,  fluctuating 
or  non  linear  part.  Velocity  variations  of  the  first  type  are  regarded  as 
signal  and  those  of  the  second  type  as  noise.  This  nomenclature  reflects 
the  hope  that  any  physical  process  of  interest  can  be  parameterized  in 
terms  of  the  large  scale  signal,  with  small  scale  noise  only  making  it 
difficult  to  measure  the  signal  veil. 

This  convention  obscures  the  true  nature  of  the  velocity  field  which 
is  that  it  has  variations  on  all  length  scales,  with  a  smooth  decrease  in 
amplitude  for  decreasing  length  scales.  There  is  no  clear  division  between 
large  and  small  scale.  Furthermore  there  is  no  clear  reason  to  associate 
large  scale  with  signal  and  small  scale  with  noise.  For  studies  involving 
the  actual  opening  and  closing  of  leads,  the  small  scale  phenomena  may  indeed 
be  the  signal  and  it  cannot  readily  be  parameterized  by  the  large  scale  motion. 

Perhaps  a  better  conceptual  model  of  the  spatial  structure  of  the  ice 
velocity  is  a  system  which  accepts  a  smooth  input  and  produces  a  discontinuous 
output,  (a  simple  example  of  such  a  system  is  the  greatest  integer  function: 

L*3  *  greatest  integer  less  than  or  equal  to  *  .) 

- The  shift  in  emphasis  from  the  earlier  model  is  this.  The  earlier  model 

viewed  the  velocity  as  the  sum  of  separate  contributions, one  smooth,  one 
discontinuous.  The  alternative  views  the  velocity  as  a  discontinuous  response 
to  a  smooth  input.  This  point  of  view  may  lead  more  naturally  to  phenome¬ 
nological  descriptions  of  the  properties  of  the  ice  pack,  through  a  comparison 
of  the  input  and  output  fields.  To  my  knowledge  this  has  not  been  attempted 
because  data  with  adequate  spatial  resolution  of  the  velocity  field  are  still 
too  scarce.  Such  a  study  would  not  explain  why  the  velocity  field  had 
certain  properties.  That  explanation  must  be  based  on  rather  deeper  understanding 
of  the  geometry  of  the  ice  pack  and  the  forces  which  act  between  floes  than  ve 
have  at  present.  Still  it  would  be  useful  to  compare  some  of  the  properties 
of  the  ice  velocity  with  properties  of  the  external  driving  forces,  i.e.,  the 
geostrophic  wind  and  the  ocean  currents. 
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5.U.£‘  ,  Errors  in  estimating  the  large  scale  deformation 

The  large  scale  deformation  estimates  of  the  previous  section  are 
subject  to  sampling  variations  since  they  depend  on  the  particular  array 
of  measurements  points.  There  have  been  several  attempts  to  quantify  the 
sampling  errors.  Hibler  et  al.,  1971*;  Nye,  1975;  and  Thorndike  and  Colony, 

1977  examined  the  departures  of  the  measured  velocities  from  the  best  fit 
linear  relationships.  These  departures,  called  residuals,  homogeneity 
variations,  or  nonlinear  fluctuations,  represent  the  variability  of  the 
actual  velocity  field  on  scales  smaller  than  the  diameter  of  the  region  of 
interest. 

From  observations  over  a  200  km  scale  taken  during  the  spring  of  1972, 
Hibler  et  al.  found  root-mean-square  velocity  residuals  of  0.06  cm  sec”  , 
compared  to  typical  linear  changes  in  velocity  (  )  over  20  kn  of 

O.lU  cm  sec**  .  Thorndike  and  Colony  used  observations  from  a  100  km  scale 
taken  during  the  spring  of  1975  and  found  rms  residuals  of  0.1*  cm  sec”  and 
typical  linear  changes  over  100  km  of  1.1  cm  sec”\  Their  summer  values  were 
somewhat  larger:  1.1  cm  sec”^  and  1.8  cm  sec  ^  for  the  nonlinear  and  linear 
contributions.  These  residuals  can  be  regarded  as  errors  with  respect  to  the 
large  scale  average  derivatives  for  the  region*  When  the  number  of  observa¬ 
tions  is  small,  the  estimated  large  scale  deformations  axe  strongly  contaminated 
by  these  errors. 

Using  what  we  know  about  the  correlation  functions  for  velocity,  we  can 
profitably  address  the  sampling  question  from  a  different  angle.  Taking  the 
line  integral  definition  for  the  large  scale  velocity  derivatives,  we  ask: 
how  many  points  around  the  perimeter  of  the  region  must  be  sampled  to  resolve 
the  integral  to  some  desired  accuracy?  Intuition  suggests  that  the  measure¬ 
ments  should  be  spaced  closely  enough  to  permit  good  interpolation  but  not 
so  closely  as  to  be  highly  redundant.  A  correlation  between  velocities  at 
neighboring  measurement  points  of  0.5  might  be  a  fair  guess;  this  would 
indicate  a  spacing  of  about  1*00  km. 

To  get  a  better  answer,  we  estimate  the  line  integral  in  equation  6 
around  a  circular  region  of  radius  r  as 
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Then  as  the  number  of  measurements  M  increases,  approaches 

Since  is  a  linear  combination  of  velocity  components,  we  can  find  its 
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variance  in  terms  of  the  variances  and  covariances  of  the  velocity  components 
at  all  M  points*  and  these  can  he  evaluated  in  terms  of  the  correlation 
functions  0|j  and  .  Let  U.  he  the  vector  involving  all  2M  velocitv 
components,  A  he  the  vector  of  coefficients,  and  ^  be  E  U.U.  .  Then 
*■  /iT U.  ££«  9  A*RA  .  This  calculation  has  been  done 

for  M  *  60  which  was  large  enough  to  resolve  .  The  dependence 

of  (£  on  the  radius  v  is  shown  in  figure  2$a  . 

The  figure  shows  that  typical  values  for  '**/>*<  are  about  1%  per  day, 
decreasing  somewhat  as  the  radius  of  the  region  increases.  Of  course  this 
curve  is  a  consequence  of  the  observed  functions  B)t  and  .  The  shape 
of  the  curve  for  radii  smaller  than  100  km  remains  speculative  until  more 
details  are  known  of  the  behavior  of  and  Ox  for  small  arguments. 

The  variance  of  the  error  in  estimating  "hu ./fe*  using  only  P)  points, 

-  tyc)1  can  calculated  in  the  same  manner  simply  by  redefining 
the  coefficients  in  the  vector  A  .  These  results  are  presented  as  a  fraction 
of  the  signal  variance  E(l^)  in  figure  2$  t  .  When  M  is  3,  as  was  the 
case  for  AIDJEX  1972,  the  1975-76  manned  AIDJEX  array,  and  LOREX  79*  and 
the  radius  was  about  100  km,  the  ratio  of  estimation  error  variance  to 
signal  variance  was  about  0.25.  For  a  radius  of  1*00  km  and  a  spacing  between 
stations  of  1*00  km,  M  *  6,  and  the  ratio  of  error  variance  to  signal  variance 
is  about  0.01.  This  confirms  our  intuition  that  a  spacing  of  1*00  km  should 
be  adequate  to  resolve  the  large  scale  average  velocity  gradients  quite  well. 

Figures  2}C  and  show  results  obtained  in  a  similar  way  for  the 
large  scale  average  derivative  'dv/ax  and  for  the  large  scale  average  diver¬ 
gence.  Accurate  estimation  of  the  large  scale  average  divergence  has  special 
importance  because  of  the  role  it  plays  in  the  theory  of  the  ice  thickness 
distribution.  If  we  want  to  keep  the  rms  error  in  divergence  below  say  20' 
of  the  rms  divergence  we  require  a  ratio  of  variances  in  figure  23J  of  0.2C  =  0.0l*. 
For  a  region  of  radius  200  kilometers  this  level  of  accuracy  can  be  achieved 
with  6  measurement  stations. 

5.1*.  ^  The  relationship  between  measured  large  scale  deformation  and  total 
_ opening  and  closing 

Suppose  a  region  of  interest  is  intersected  by  a  number  of  cracks,  each 
of  which  is  opening  or  closing  at  some  rate.  If  the  rates  were  known,  the 
total  rates  of  opening  and  closing  for  the  region  could  be  found.  In  practice 
we  are  not  able  to  measure  the  motion  at  each  crack,  but  only  the  motion  of  a 
few  points  in  the  region.  How  are  we  to  use  these  few  measurements  to  esti¬ 
mate  the  opening  and  closing? 
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In  an  analogous  situation  in  one  dimension,  ve  have  k  cracks  each  vlth 
opening  or  closing  rate  Ufc  i* !,•••,  k  ♦  We  can  imagine  the  velocity  k  to 
he  a  random  function  of  having  random  discontinuities  at  random  points 
and  being  constant  between  the  points.  Suppose  ve  have  measured  only  the 
motion  at  the  end  points  of  the  region,  an  interval  of  length  L  •  Then 

A  *  4*VoJl  • •) 

'  J  *.») 

Q  4*V>I  s  (*it  *) 

*  4*  1 

k 

s  A  +  C  •  X  “>■ 

and  the  problem  is  to  estimate  A  and  C  given  U  .  At  first  sight  the 
situation  seems  hopeless.  U.  clearly  contains  information  about  the  net  open¬ 
ing  or  closing  but  not  about  the  total  opening  and  closing.  However,  knowledge 
about  the  random  variables  If  u,,Wt/ •  can  be  used  to  make  probabalistic 

statements  about  the  opening  and  closing.  Suppose  for  instance  k  has  the 
Poisson  distribution  with  parameter  X  .  This  means  that  is  the  average 

number  of  cracks  in  an  interval  of  length  L  .  Then  the  orobability  of 
finding  k  cracks  in  a  random  interval  of  length  L.  is 

rfM.L)  •  • 

If  >L  -  0.1,  for  instance,  the  probability  of  getting  k  *  0  is  0.9051, 

*  0.0905,  and  *  n.00^5.  Therefore,  with  high  probability, 

there  is  either  no  crack  or  Just  one  in  the  interval  and  the  observed  value 

of  the  velocity  difference  ll  would  itself  indicate  the  total  opening  and 

total  closing.  Of  course  for  larger  it  becomes  more  likely  that  several 

cracks  intersect  the  test  interval,  in  which  case  the  observed  U.  cannot 

separate-  the  opening  from  the  closing. 

A  similar  approach  for  the  two  dimensional  problem  is  to  imagine  the  ice 

pack  to  be  crisscrossed  by  a  family  of  random  cracks,  defined  by  the  random 
straight  lines 

%  cmBi  ♦  ^  *  rc 

where  each  has  the  uniform  distribution  on  (b,27)  and  the  normal  distances 


1 

(7) 
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from  the  origin  to  each  line  form  a  sequence  of  Poisson  points.  Imagine 
that  associated  vith  each  crack  is  a  velocity  discontinuity  having 
the  Gaussian  distribution.  This  describes  a  random  vector  field  having 
structure  akin  to  the  velocity  field  of  sea  ice.  The  observations  of  Hall 
and  Rothrock  can  be  used  to  estimate  the  parameters  in  the  Poisson  and 
Gaussian  distributions.  These  random  fields  have  the  following  properties. 
They  consist  of  discrete  rigid  floes.  The  floes  have  a  distribution  of  sizes 
determined  by  the  Poisson  field  of  lines.  The  Poisson  lines  are  isotropic 
and  homogeneous.  The  velocity  difference  between  any  two  points  is  the 
vector  sum  of  the  velocity  discontinuities  encountered  getting  from  one  point 
to  the  other. 

At  each  crack  the  opening  or  closing  is  determined  by  the  projection 
of  the  velocity  difference  vector  onto  the  normal  to  the  crack. 


opening  *  max  (  o , 

c&j  +  v;  ) 

closing  ■  min  (  o, 

<*L’vv  ) 

Thus  it  is  a  simple  matter  to  evaluate  the  total  opening  and  closing  for  any 
realization  of  the  random  field. 

Of  course  one  can  also  imagine  measuring  the  velocity  at  a  few  points 
and  constructing  the  /.-strain  rate  invariants  from  the  observed  velocities. 
In  this  way  one  can  test  for  a  relationship  between  the  -strain  rate 
invariants  calculated  from  a  few  sampled  velocities  and  the  total  opening 
and  total  closing  found  by  tallying  up  the  activity  at  every  crack. 

In  an  attempt  to  carry  out  this  program,  parameter  values  were  taken 
to  be  A  *0  0? km”1,  A  *  ^  (arbitrary  units)2.  An  L  *  100  kr. 

triangle  was  used  to  simulate  the  sampling  procedure  used  during  AIDJEX. 

A  large  number  of  random  fields  were  generated. 

Each  realization  of  the  random  field  is  defined  by  the  sequences  of 
random  numbers  k ,  &i  ,  tt;  ,  V,;  #  f;  f or  .  Here  the  0L  are  drawn 

independently  from  the  uniform  distribution  on  .  The  values  for 

form  a  Poisson  process  with  parameter  ^  .  This  is  achieved  by  drawing 
the  increments  A*  1^, |  -T^  independently  from  the  exponential  distribution 

with  density  .  The  process  is  terminated  as  soon  as  exceeds 

100  kilometers  since  none  of  the  subsequent  lines  would  intersect  the  100 
kilometer  region.  Finally  the  and  V;  were  drawn  independently  from 
the  normal  distribution  with  zero  mean  and  unit  variance. 
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For  each  field  the  velocity  was  measured  at  the  three  specified  points 
and  the  -strain  rate  invariants  and  &  were  evaluated. 

Also  for  each  field  the  total  opening  and  closing  vere  evaluated  using  eq.  7  • 

These  vere  normalized  by  |£]  and  displayed  versus  9  in  figure  li ,  For  each 
realization  of  this  random  field  two  points  are  plotted, (  B,  A/\Fi)  and 

)  •  From  the  plot  it  is  clear  that  there  are  not  unique  values 

of  A/iE}  end  C/iEI  corresponding  to  a  given  B  .  Instead  there  is  a 
distribution  of  values  for  A/ F|  (and  for^EI  ),  and  this  distribution 
changes  with  B  .  The  distribution}  sketched  in  figure  2$ t  sure  broad  in  the 
sense  that  probable  departures  from  the  mean  are  at  least  as  large  as  the  mean. 

The  interpretation  of  this  exercise  is  that  the  total  opening  and  closing 
are  only  weakly  determined  by  the  t.  -strain  rate  invariants  based  on  three 
measurements  100  kilometers  apart.  Had  the  relationship  been  a  strong  one 
the  distributions  in  figure  25  would  have  been  narrower,  or  to  say  it  differ¬ 
ently,  the  points  in  figure  2f  would  have  clustered  more  closely  around 

curves  like  those  in  figure  12  .  It  may  still  be  useful  to  imagine  these 

smooth  curves  but  only  with  the  recognition  that  the  actual  opening  and 

closing  scatter  widely  about  the  imagined  curve: 

*  Wa(  B)  +  random  error 

C/  ,  ♦  random  error 

/\£)  • 

where  the  random  terms  and  the  6 - dependent  terms  make  roughly  equal  contri¬ 
butions  to  the  total  opening  and  closing. 

The  total  opening  and  closing  are  an  essential  part  of  the  theory  of 
the  ice  thickness  distribution.  The  above  results  suggest  that  it  will  be 
difficult  to  estimate  the  opting  and  closing  accurately.  Perhaps  the  best 
that  can  be  hoped  is  to  drive  the  thickness  distribution  calculations  with 
opening  and  closing  time  series  which  have  the  right  statistical  properties 
even  if  they  may  have  large  errors  on  a  day  to  day  basis.  The  statistical 
properties  can  be  inferred  from  figure  2V  .  Of  course  this  figure  is  the 

consequence  of  a  particular  conceptual  model  motivated  by  a  limited  data 

set— the  SEASAT  SAP  data  from  early  October  1978.  More  data  need  to  be 
collected  and  studied  before  these  ideas  can  be  extended  to  other  times  and 
places  with  confidence. 
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ionship  'between  kinematics  and  stress 
Kinematic  data  can  be  used  for  studying  the  relationship  between 
stress  and  strain  for  sea  ice.  Suppose  for  instance  that  the  state  of  stress 
tensor  (T  suitably  averaged  over  some  region ,  is  related  to  kinematic  quan¬ 
tity*  c  and  the  ice  geometry  5  by  some  expression  of  the  form 

<T  *  FU.s)  . 

Although  <r  cannot  be  measured  directly,  its  divergence  V-T  can  be  inferred 
indirectly  from  the  observed  momentum  balance  when  all  the  other  terms  in 
that  balance  are  known.  If  a  particular  function  F  is  hypothesised,  it 
can  be  tested  by  evaluating  both  sides  of  the  equation  *  V-  r(t,s)  t 

the  left  hand  side  as  a  residual  from  the  momentum  balance  and  the  right 
hand  side  in  terms  of  observed  kinematic  and  ice  geometry  quantities.  This 
was  one  of  the  objectives  of  AIDJEX  (Maykut  et  al.,  1972).  In  practice 
the  test  is  difficult  to  make  since  neither  V  (T  nor  can  be  deter- 

mened  very  accurately  from  observations.  See  Bothrock  et  al.,  1980. 

Under  special  circumstances,  meaningful  tests  may  be  possible,  however. 
When  the  ice  deformation  is  strongly  divergent,  ice  floes  tend  to  move  apart. 
With  no  floe-to-floe  contacts  there  can  be  no  V-  0*  forces.  The  left  hand 
side  should  differ  from  zero  only  by  measurement  error.  These  errors  are 
probably  small  enough  to  permit  a  useful  test. 

Another  special  situation  of  interest  is  when  the  ice  is  being  forced 
up  against  the  coasts  by  the  wind.  As  the  ice  moves  toward  the  coast  it 
converges,  becomes  stronger,  and  eventually  becomes  strong  enough  to  resist 
further  deformation.  If  the  on  shore  winds  persist,  a  zone  of  motionless 
ice  can  widen  to  several  hundred  kilometers  (Pritchard,  1977).  In  this  situa¬ 
tion,  intuition  suggests  that  the  ^*0"  vector  should  point  off  shore  and 
should  increase  in  time  until  it  balances  the  on  shore  wind  stress.  With 
simultaneous  ice  trajectories  at  100,  200,  300  and  1*00  km  from  the  coast 
it  should  be  possible  to  observe  the  amount  of  convergence  required  to  pro¬ 
duce  the  required  resistance  to  the  wind. 

At  greater  distances  from  th  coasts,  there  is  evidence  that  the  ice 
stresses  embodied  in  V*<T  are  usually  small.  Little  success  is  anticipate 
in  trying  to  observe  them.  Any  stress-strain  law  which  provides  adequate 
resistance  to  deformation  near  the  shore  is  probably  adequate  for  full  basin 


dynamic  modelling.  The  observations  may  never  be  adequate  to  discriminate 
very  selectively  betveen  candidate  stress-strain  lavs. 

5.5  Discussion 

I  have  tried  in  this  paper  to  bring  together  some  of  what  is  known 
about  the  motion  of  sea  ice.  The  emphasis  has  been  on  the  departures  of 
the  true  motion  from  the  long  term  mean  circulation.  This  is  a  compari- 
tively  new  topic,  the  investigation  of  which  has  only  recently  been  made 
possible, mainly  by  the  developments  of  satellite  positioning  techniques. 
Perhpas  because  of  its  newness,  the  field  lacks  a  clear  agenda  of  questions 
to  be  answered.  Instead  we  are  still  trying  to  characterize  the  motion,  to 
determine  the  magnitude  of  the  motion  on  different  scales  and  to  identify 
motions  which  have  some  intrinsic  interest  or  are  related  to  other  phenomena 
of  interest.  It  is  not  surprising  that  contributors  to  the  field  have  dif- 
ferred  in  their  approaches  to  observing  and  interpreting  the  motion. 

The  time  and  space  corelations  have  been  used  extensively  here  for 
several  reasons.  First  they  give  a  compact  description  of  the  motion.  In 
principle,  the  correlation  function  involves  one  time  and  two  space  variables, 
but  in  practice  a  great  deal  of  information  can  be  extracted  from  three 
functions , each  of  a  single  variable-  Gufr)  and  .  Second, 

these  three  functions  are  accessible  to  observation.  Third,  prope  ' ies  of 
nearly  all  kinematic  quantities  can  be  deduced  from  these  functions, as  illus¬ 
trated,  ad  nauseura,  in  the  text.  Fourth,  these  functions  form  the  rigorous 
basis  for  answering  questions  related  to  interpolation,  prediction,  and 
experiment  design. 

Many  of  the  results  presented  here  are  based  on  sample  autocorrelation 
functions  deduced  from  limited  data.  The  data  available  for  esti¬ 

mating  correlations  at  small  space  lags  (less  than  100  km)  are  meagre  indeed. 
More  work  along  the  lines  of  Hall  and  Rothrock  would  help  to  resolve  this 
part  of  the  correlation  function.  The  behavior  of  the  correlation  functions 
G  (<r)  o.*}  G|_(i-)  in  the  limit  of  small  r  is  an  important  property  of 

the  motion  related  to  the  granular  nature  of  the  ice  pack. 

As  mentioned  in  the  text, care  must  be  exercised  in  choosing  correlation 
functions  or  contradictions  (negative  variances)  can  occur.  In  fact  this 
has  happened  in  some  of  the  calculations  done  using  the  correlations  tabulated 
in  Tabled  with  linear  interpolation  to  intermediate  distances.  This  means 
that  the  piecewise  linear  function  defined  in  the  table  is  not  positive 
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VI.  Spatial  correlation  functions  for  sea  ice  v 
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Figure  5.  (a)  Latitude  of  ice  station  Caribou,  versus  time,  2  Sept-16  Nov  1975. 

NavSat  data,  sampling  interval,  2.5  days. 

(b)  Latitude  of  Ice  station  Dig  Bear,  versus  time,  2-15  Sept  1975. 
Navoat  data, sampling  interval,  to  hour;;. 

(r)  Arbitrary  y-ooord  1  nake  of  Dig  Hear,  versus-  time,  ?-1  fept  1 97 '> • 
Ai'oush i  r  trai'king.  data,  sampling  Interval  15  minute:;. 

(d)  Arbitrary  x-coordinute  of  Dig  Buur,  versus  time,  ?.  Cept  1975. 
Acoustic  tracking  data,  sampling  interval  one  minute. 

Only  in  (d)  is  the  measurement  error  visible  (±3  m). 
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Figure  11.  The  complex  time  correlation  function  for  ice  station  Caribou,  1975-76. 

Trie  velocity  variance  is  1^5  cm  2  r,ec“  *. 

Figure  Ik1.  The  time  correlation  function  for  the  u  and  v  components  of  velocity 

from  28  grid  points  in  the  central  part  of  the  Arctic  Basin,  from  drifting 
buoy  data  collected  during  1979.  The  cross  correlation  between  u.  and 
v  in  negligible.  Variance  of  U.  =  P3  cm*  sec-1  ;  variance  of  v  *  ?2  cm1  see- 
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Figure  13.  (a)  The  real  and  imaginary  parts  of  the  weight  a*,  to  be  used  for  inter¬ 

polating  to  the  mid  point  of  a  time  interval  of  duration  X  .  The  ratio 
of  the  measurement  error  to  the  standard  deviation  of  velocity  is  cr/£ 
(b)  The  variance  of  the  interpolation  error  expressed  as  a  fraction  of 
the  variance  of  the  ice  velocity.  Data  from  figure  11  were  used  hero. 
Logarithmic  time  scale. 
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Figure  lU.  (a)  The  weight  ex'  for  prediction  over  time  t  . 

(b)  The  prediction  error  variance  versus  prediction  time.  See  Figure  13 
also. 


Figure  15.  A  schematic  representation  of  the  X  and  £  components  of  velocity  at 
the  points  (  o,t  ),  (r,o)t  (  e,r).  Isotropy  implies  EaC  *  E  Wl( 

£Wi  '  Eat  »  Mid  •  EJf  t  but  does  not  imply  equality  between 
tac  end  Eae  or  Elf  and  £  11 
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Figure  16.  Observed  correlations  between  velocity  components  &'•'  functions  of  di stance. 


Figure  17.  The  variance  of  the  interpolation  error  as  a  function  of  the  grid  spacing, 
and  the  measurement  error,  o*1  . 


Figure  18.  The  strain  rate  invariant  0  indicates  whether  the  motion  is  dominated 
by  divergence,  shear,  or  convergence.  The  total  opening  and  closing 
are  sometimes  treated  as  functions  of  P. 
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Figure  19- 
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The  variance  of  velocity  increments  versus  interval  length,  from  SeaSat 
synthetic  aperture  radar  data,  3-5  October  1978,  Beaufort  Sea. 


Figure  20.  Histograms  of  the  strain  rate  invariant  0  .  Solid  line  is  for  drifting 

buoy  data  from  1979,  grid  spacing  about  500  km  (Colony  and  Thorndike,  1901) 
Dots  are  for  AIDJEX  manned  camp  data,  spacing  100  km  (Hye,  1976). 
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Figure  21.  Strain  ellipses  for  the  AID.7EX  buoy  array,  800  km  diameter  (a);  and  for 
the  AIDJEX  manned  camp  array,  200  km  diameter  (b).  The  ellipses  shov 
the  deformation  of  a  circular  region  on  1  May  1975,  to  the  date  indicated. 
The  principal  axes  of  the  ellipse  are  the  principal  strains.  The  angle 
from  the  horizontal  broken  line  to  the  major  axis  is  the  princinal 
direction.  The  rigid  body  rotation  is  indicated  by  the  arc  f rom  rn  to  «. 
Data  from  the  Beaufort  Sea,  roughly  a  500  km  radius  about  7i»°N,  ll*5c'W. 
(From  Thorndike  and  Colony,  1980.) 


Figure  22.  Strain  ellipses  from  1979  (solid  line)  and  I960  (broken  line)  driftinr 
buoy  data,  showing  the  year  long  deformation  of  an  initial  circle 
(drawn  over  Greenland). 
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Figure  23.  (a)  Typical  values  for  large  scale  average  strain  rate  quantities  are 

functions  of  the  radius  of  the  region  over  wnicli  the  average  is  taken, 
(b)  The  estimation  error  for  various  average  strain  rate  quantities 
depends  on  the  radius  of  the  region  and  on  the  number  of  measurements 
made  around  the  perimeter. 
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Figure  Tin*  i.oljil  opening  and  closing  versus  3  .  10  realizations  of  the  random 

field  wii.li  I'ol-'.r.oti  parameter  AL  *  B.  The  velocity  discontinue  ties 
were  drawn  independently  from  the  unit  normul  distribution. 
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Figure  25.  Distributions  of  the  total  opening  for  several  ranges  of  9  .  Data  take 
from  figure  2 1*. 


